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I. INTRODUCTION 

The possibility suggested by Weinberg fil and pio- 
neered by Ray, Ordonez and Van Kolck [1115 of making 
model independent predictions for NN scattering using 
Effective Field Theory (EFT) methods and, more specif- 
ically, Chiral Perturbation Theory (ChPT) has triggered 
a lot of activity in recent years (for a review see e.g. 
Ref. p|)- In addition to the previous works, most subse- 
quent calculations dealing with the specific consequences 
of ChPT have been focused in making predictions for 
NN scattering phase-shifts and deuteron properties based 
in thegenuine Two Pion Exchange (TPE) chiral poten- 
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HO, HH HI, HI, HI, HI, HI, Ha although some incipient 
work has also recently been started imp lementing Three 
Pion Exchange effects 0, 0, 0, El El] . In a given 
partial wave (coupled) channel with good total angular 
momentum the reduced NN potential (U(r) = MV(r)) in 
configuration space can schematically be written in local 
and energy independent form |7l llllll2l | (see Refs. |2j.y.|_l( 
for an energy dependent representation) for any distances 
larger than a finite short distance radial regulator, r c , 
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-j^Wm!ho(mr, g, a, c 3 , c 4 ) 



(1) 



where W(x) are known dimensionless functions which are 
everywhere finite except for the origin where they ex- 
hibit power law divergencies, which demand the use of 
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some regularization. In writting the previous expression 
we have disregarded distributional contact terms (deltas 
and derivatives of deltas) which strengths are scheme de- 
pendent but do not contribute for r > r c > 0. The po- 
tential is completely specified by the pion mass, m, the 
pion weak decay constant, /, the nucleon mass M, the 
axial coupling constant g, the Goldberger-Treimann dis- 
crepancy dis and three additional low energy constants 
ci = C\M, C3 = C3M and C4 = C4M which can directly 
be deduced from the analysis of low energy ttN scatter- 
ing within ChPT |M El El E| . Given this informa- 
tion one can then solve the single or coupled channel 
Schrodinger equation imposing a regularity condition of 
the wave function at the origin for each separate chan- 
nel. In this paper we will work under the assumption 
that the long range pieces of the potential should be it- 
erated to all orders, but some perturbative analysis will 
also be done. Our motivation is to describe long range 
correlations between observables in the NN problem in 
a model independent way. As we will see below there is 
still some additional physical information required in the 
form of either counter-terms or short distance boundary 
conditions to make the problem well posed if one indeed 
wants to go to remove the regularization. They depend 
exclusively on the short distance behaviour of the poten- 
tial through phases of the wave function. The number 
of phases depends crucially on the repulsive or attractive 
character of the potential. In this sense the power count- 
ing for the short distance interactions cannot be regarded 
as independent on the power counting of the singular 
chiral potentials. Non-perturbatively this materializes, 
after renormalization, in non-integer power counting for 
physical observables. This imposes severe limitations on 
the admissible structure of counterterms and the corre- 
sponding renormalization conditions of the quantum me- 
chanical problem. Right away we hast to emphasize that 
our approach is not the conventional EFT one of allow- 
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ing all possible short distance counterterms allowed by 
the symmetry, and to a certain extent our viewpoints 
are admittedly heterodox within the conventional EFT 
framework. However, based on the physical requirement 
of having small wave functions in the short range un- 
known region the basic and orthodox quantum mechan- 
ical requirements of completeness and orthoghonality of 
states are deduced, providing a justification for the ad- 
ditional restrictions. We remind here the series of works 
by Phillips and Cohen ( see also R e - H3)> where 

restrictions on zero range interactions were deduced for 
non-singular potentials based on the Wigner causality 
conditions. Here we extend their results also to the sin- 
gular NN interactions of Eq. . 

The theorem underlying the EFT developments is that 
if chiral symmetry is spontaneously broken down in QCD, 
the true NN potential at long distances is embedded in 
the parameter envelope of the general chiral NN poten- 
tial, Eq. ifTJl. and the chiral expansion provides a re- 
liable hierarchy at those long distances. The hope is 
that compatible and perhaps accurate determinations of 
both ttN and NN low energy data, bound states and 
resonances can be achieved with the same sets of pa- 
rameters. The problem is that, in order to make truly 
model independent predictions, short distance ambigui- 
ties should be under control and their size smaller than 
the experimental data uncertainties used as input of the 
calculation. Only then can the renormalization program 
be carried out satisfactorily as it was done in the OPE 
case |33, Ha EJ EH UM > although, as recognized by 
Nogga, Timmermans and Van Kolck, this may be done 
at the expense of modifying the power counting 42] of 
the counterterms in favor of renormalizability (see also 
Ref. [41j for a complementary formulation in terms of 
boundary conditions). The present work analyses this 
problem extending our previous OPE renormalized cal- 
culations to NNLO TPE and its implication in the values 
of the chiral constants. 

The determination of the chiral constants c\ , c 3 and C4 
(in units of GeV" 1 from now on) from nN scattering has 
been undertaken in several works and shows significant 
systematic discrepancies depending on the details of the 
analysis. In Heavy Baryon ChPT for low energy irN scat- 
tering 30] the values a = -1.23±0.16, c 3 = -5.94±0.09 
and C4 = 3.47 ± 0.05 were deduced with a sigma term of 
<t(0) = 70MeV. In Ref. |3lJ the analysis of low energy 
ttN scattering inside the Mandelstam triangle yields c\ = 
-0.81±0.15, c 3 = -4.69± 1.34 and c 4 = 3.40±0.04 with, 
however, a bit too low sigma term tr(0) = 40MeV as com- 
pared to ChPT. Unitarization methods reproducing the 
phase-shifts [33, HH from threshold to the A resonance 
region conclude a = -0.43±0.04, c 3 = -3.10±0.05 and 
c 4 = 1.51 ±0.04. 

The values of the chiral constants ci, C3 and C4 also 
depend on regularization details of the NN chiral in- 
teraction. The irN values from Ref. [3lJ were taken in 
the nucleon-nucleon NNLO calculation of Ref. 01 with 
sharp and gaussian cut-offs A = 0.6 — 0.8GeV in mo- 



mentum space, and momentum dependent counter-terms 
were supplemented and determined from a fit to the NN 
data base based on the Partial Wave Analysis (PWA) 
of Ref. [HEil. Likewise, Ref. constructs a NNLO 
chiral potential where channel dependent gaussian mo- 
mentum space cut-offs in the range A = 0.4 — 0.5GeV 
were used to fit the NN database HH. The N 3 LO exten- 
sion of this work |28| uses only one common cut-off and 
fixing Ci = —0.81 produces c 3 = —3.20 and C4 = 5.40. In 
Ref. [ll| the NNLO calculation was done in configuration 
space with a short distance cut-off at r — 1.4fm where an 
energy and channel dependent boundary condition was 
imposed and the fixed value c\ — —0.76 ± 0.07 was used 
to make a PWA to pp data yielding c 3 = —5.08 ± 0.24 
and C4 = 4.70 ± 0.70. An update of this calculation also 
including np data [2l| generates c 3 = —4.78 ± 0.10 and 
c 4 = 3.96 ± 0.22. The calculations of Ref. ^ im- 
prove the cut-off dependence of the potential in momen- 
tum space by using spectral regularization and take again 
gaussian cut-offs and fix ci = —0.81 yielding after fitting 
the counter-terms to the NN PWA Oj , the values 
c 3 = —3.40 and C4 = 3.40. The extension of this work to 
N 3 LO has been done in Ref. [2^] keeping the same values 
for c 3 and C4 and readjusting the counter-terms. 

In a renormalized theory results should be insensitive 
to the auxiliary regularization method if the regulator is 
removed at the end. If a fit to the database proves suc- 
cessful, then the resulting parameters should be cut-off 
independent or at least the systematic uncertainty in- 
duced by the regularization should be smaller than the 
statistical errors induced by experimental data. Oth- 
erwise, the cut-off becomes a physical parameter. The 
first indication that finite cut-offs effects are sizeable in 
present calculations has to do with the variety of values 
that have been used in the literature for the low energy 
constants c±, c 3 and ca to adjust N N p artial waves and 
deuteron properties [13, IH II |H HI H| (see also the 
comment in Ref.[24|]). Obviously, we do not expect the 
values of the c's to agree exactly, but the discrepancies 
should be at the level of the difference in the approxima- 
tion . Since the data base is the same but the regular- 
ization schemes are different, one unavoidably suspects 
that these determinations of the LECS may perhaps be 
regularization and hence cut-off dependent. 

To get a proper perspective on the issue of renormaliza- 
tion let us consider the size of the contributions of the chi- 
ral potential in configuration space at different distances. 
For instance, at r = 1.4fm in the 1 Sq channel each order 
in the expansion is about an order of magnitude smaller 
than the preceding one. At short distances, however, the 



For instance, pion loops at NLO modify the contribution to C3 
by ~ 3g^m 2 /(647r/ 2 ) ~ 0.4/GeV. This contribution must be 
tak en into account when comparing numbers between lilt |2 ll , 
l.'ill and the present approac h. Only the extractions using the 
N3LO NN pot ential in Ref. |28l| are made at the same order as 
those from l.'ill . 
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situation is exactly the opposite, higher orders dominate 
over the lower orders. In the previous example of the 1 Sq 
channel, LO and NLO become comparable at r ~ 0.9fm, 
and NLO and NNLO become comparable at distances 
which value r ~ 0.1 — 0.4fm depends strongly on the par- 
ticular choice of low energy constants C3 and C4. Actually, 
a general feature of the chiral NN potentials at NNLO has 
to do with their short distance behavior; they develop 
an attractive Van der Waals singularity U ~ —MC^/r 6 
similar to the one found for neutral atomic systems. In 
such a situation, the standard regularity condition at the 
origin only specifies the wave function uniquely if the po- 
tential is repulsive, but some additional information is 
required if the potential is attractive 0| (for a compre- 
hensive review in the one channel case see e.g. Ref. [47j.L 
Within the EFT framework the problem has been revis- 
ited in Ref. 0]. The net result is that the regularity 
condition at the origin tames the singularity |4l| and, in 
fact, more singular potentials become less important at 
low energies. 

In this work we reanalyze the NN chiral potential in- 
cluding TPE potential at NNLO. We carry out the anal- 
ysis entirely in coordinate space following the ideas de- 
veloped in our previous work |4l| for the OPE potential. 
In configuration space the (renormalized) potential is fi- 
nite except at the origin, a point which should carefully 
be handled, requiring a delicate numerical limiting pro- 
cedure. For a singular potential in coordinate space, the 
corresponding potential in momentum space is not finite 
unless a short distance cut-off or a subtraction proce- 
dure at the level of the potential is implemented, hence 
modifying the potential everywhere and not just at high 
energies. This results generally in two cut-offs: one for 
the irreducible two point function and another for the 
Lippmann-Schwinger iteration [Tol. [T9II20II2I I2II2I l29j . 
The short distance coordinate space cut-off is unique and 
common both to the potential and the scattering solu- 
tion. 

Unlike previous works on the TPE potential we try to 
remove the cut-off completely taking the consequences 
seriously. This does not mean that finite cut-off calcula- 
tions are necessarily incorrect or not entitled to describe 
all or some part of the data, but there are also good rea- 
sons for removing the cut-off and looking at the physical 
consequences. In the first place the limit exists in strict 
mathematical sense under well defined conditions, as the 
analysis below shows. This is a non-trivial fact, because 
calculations done in momentum space can only address 
this question numerically by adding counterterms sug- 
gested by an a priori power counting on the short dis- 
tance potential. As shown in Ref. jij this does not al- 
ways work, and calculations may require some trial and 
error. Secondly, this is the only way we know how to 
get rid of short distance ambiguities, and thus to make 
truly model independent calculations. Third, the study 
of peripheral waves has proven to be successful by using 
perturbative renormalized amplitudes corresponding to 
irreducible TPE and iterated OPE where the cut-off has 



been removed in the intermediate state I! l2a. Periph- 
eral waves mainly probe large distances in the Born ap- 
proximation but they also see some of the short distance 
interaction due to re-scattering effects. Fourth, the ad- 
vantage of renormalization is that one should obtain the 
same results provided one uses as input the same physical 
information, regardless whether the calculation is done in 
coordinate or momentum space, and also regardless on 
the particular regularization. Finally, a reliable estimate 
on the errors and convergence rate of the chiral expansion 
can be done, without any spurious cut-off contamination. 
In principle, the higher order in the chiral expansion the 
better, provided there is perfect errorless data to fit the 
increasing number of low energy constants appearing at 
any order. However, the chiral expansion may reach a 
limited predictive power because of finite experimental 
accuracy in the low energy constants used as input. The 
output inherits a propagated error which may eventu- 
ally become larger than the experimental uncertainty 2 . 
Finite cut-off uncertainties are not a substitute for propa- 
gating input experimental errors to the predictions of the 
theory, and can be regarded at best as a lower bound on 
systematic errors. In this paper we regard this possible 
cut-off dependence as purely numerical inaccuracies of 
the calculation, and not as a measure of the uncertainty 
in the predictions of the theory, so we make any effort to 
minimize these cut-off induced systematic errors. 

In the process of eliminating the cut-off we find some 
surprises, and effects not explored up to now become 
manifest. Even for low energy scattering parameters 
and deuteron properties, where the description should 
be more reliable and robust we find systematic discrep- 
ancies in our calculation with values quoted in the liter- 
ature and which we conclusively identify as finite cut-off 
effects. This might provide a natural explanation why 
calculations with different cut-off methods fitting the NN 
phase shifts [43, 144L |45| obtain different results for the 
chiral constants c\, C3 and C4 or why different values of 
the constants yield good fits to the data. According to 
our study, for the lowest phases the reason can partly be 
related to the dominance of short distance Van der Waals 
singularities for a system with unnaturally large scatter- 
ing lengths or a weakly bound state as it is the case for 
the and 3 Si — 3 D\ channels. In some cases, they are 
as large as a 30% effect like in the effective range of the 
triplet 3 Si channel. The size of the effect depends on 
the value of the low energy irN constants ci, C3 and C4. 
Given the significant sensitivity of low energy NN prop- 
erties and deuteron properties on these low energy ttN 
constants we try to make a fit to some low energy prop- 
erties which uncertainties are reliably known and where 
we expect the chiral theory to be most reliable. At this 



2 This issue has been illustrated in Refs. I -111 15(1 l5ll for the case 
of TV-it scattering at two loops, and will become clear in NN scat- 
tering below. 
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point we depart from the standard large scale fits to all 
phase shifts or partial wave analysis where the low en- 
ergy threshold parameters are determined a posteriori. 
The assignment of statistical errors on the fitting param- 
eters ci , C3 and C4 is often not addressed (see however 
Refs. [ill l2ll|) b ecause the NN data bases used to fit the 
phase shifts |43l El Efij are treated as errorless. We also 
try to improve on this point within our framework. 

The paper is organized as follows. In Sect. [H] we dis- 
cuss the basic assumption of the smallness of the wave 
function in the short range unknown region and its conse- 
quences. We also analyze the constraints based on causal- 
ity and analyticity of the S— matrix. In Sect. Hill we in- 
troduce the classification of boundary conditions which 
will be used along the paper to effectively renormalize the 
amplitudes both in the one-channel as well as the cou- 
pled channel case. We will also review the orthogonality 
constraints for singular potentials already used in our 
previous work for the OPE potential. Sect. ITvl deals 
with the description of the singlet channel. From the 
superposition principle of boundary conditions we show 
how a universal form of a low energy theorem for the 
threshold parameters as well as for the phase-shift arises. 
In Sect. 0we discuss the interesting triplet 3 Si — 3 D\ 
channel both for the deuteron bound state as well as the 
corresponding scattering states, where full use of the or- 
thogonality constraints as well as the superposition prin- 
ciple of boundary conditions generates interesting analyt- 
ical relations connecting deuteron and scattering proper- 
ties. In Sect. IVII a careful discussion of errors for our 
cut-off independent results is carried out. Also, a de- 
termination of the chiral constants based on low energy 
data and deuteron properties is made. In Sect. IVIII we 
present a simplified study on the significance of the chi- 
ral Van der Waals forces and the striking similarities with 
the full calculations for the s-waves. In Sect. IVIIII we 
show some puzzling results for the NLO calculation in 
the deuteron channel. We also comment the relation to 
finite cut-off calculations and the conflict between Wein- 
berg counting and non-perturbative renormalization at 
NLO. We also outline possible solutions to this problem. 
In Sect. IIXI we analyze our results on the light of long 
distance perturbation theory, reinforcing the usefulness 
of non-perturbative renormalization due to an undesir- 
able proliferation of counterterms. Finally, in Sect. Ixl we 
summarize our conclusions. 

For numerical calculations we take / w = 92.4MeV, 
m = 138.03MeV, M = M p M n /(M p + M n ) = 
938.918MeV, g„ NN = 13.083 in the OPE piece to account 
for the Goldberger-Treimann discrepancy according to 
the Nijmegen phase shift analysis NN scattering and 
c]A = 1-26 in the TPE piece of the potential. The val- 
ues of the coefficients c\ , C3 and C4 used along this paper 
are listed in Table |U for completeness. The potentials in 
configuration space used in this paper are exactly those 
provided in Ref. @, [H d but disregarding relativistic 
corrections, M/E — > 1. 



II. SHORT DISTANCE INSENSITIVITY 
CONDITIONS AND RENORMALIZATION 

In this section we ellaborate on the essential role played 
by standard quantum mechanical orthogonality and com- 
pleteness properties of the wave functions in the rest 
of this paper. As we have already mentioned, our ap- 
proach is unconventional from an EFT perspective, and 
at present it is unclear whether such properties have an 
EFT justification. At the same time one should say that 
many self-denominated EFT calculations do indeed nor- 
malize deuteron wave functions to unity and use energy 
independent regulators from which orthogonality rela- 
tions follow automatically. 



A. The inner and the outer regions 

Similar to EFT, our basic assumption is that low en- 
ergy physics should not depend on short distance fine 
details. This rather general principle can be made into a 
precise quantitative statement in practice for a quantum 
mechanical system. For the sake of clarity let us consider 
the singlet 1 Sq channel for positive energies. If we as- 
sume a short distance regulator r c , above which our long 
distance (local) potential acts, the reduced Schrodinger 
equation in the outer region reads 

- u 'Lh( r ) + Uh{r)u k ^(r) = k 2 Uk,h(r) , r > r c (2) 

where the label L stands for long. Asymptotically be- 
haves as 

u k ,h(r) -> Asm(kr + S(k,r c )) (3) 

where A is an arbitrary normalization constant and the 
dependence on the short distance regulator r c has been 
explicitly highlighted. In the inner region, the dynamics 
is unknown but we also expect it to be irrelevant provided 
kr c -C 1, i.e. if we assume the corresponding wavelength 
to be larger than the short distance scale. The poten- 
tial can be deduced from perturbation theory in the full 
amplitude, 

U{x) = C S(x) + C 2 {V 2 ,<5(f)} + • • • + U L {x) (4) 

where Ul(x) corresponds to the expansion in Eq. □ 
ITU IT^ . The distributional contact terms are regular- 
ization scheme dependent and correspond to polynomial 
terms in momentum space. Obviously, they do not con- 
tribute to the region r > r c for r c > 0. The very nature 
of such a calculation already implies that Co, C%, etc. are 
perturbative corrections to the short range physics, but 
do not include possible non-perturbative effects. As will 
become clear below (Sect. ITTTt . finiteness of the physical 
phase shitf in the limit r c — ► 0, implies a highly non- 
perturbative reinterpretation of the short range terms, 
even if the long range pieces are computed perturbatively. 
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Following |3 it is useful to use a nonlocal and en- 
ergy independent potential to describe the short distance 
dynamics 

- u ks( r )+ / Us(r,r')uk,s(r') = k 2 u k ,s{r) , 
Jo 

r < r c (5) 

where the label S stands for short. This holds provided 
we are below any inelastic channel such as ttNN. Above 
such threshold, any open channel should be included ex- 
plicitly. The nonlocal short distance potential Us(r,r') 
encodes, in particular, contact terms (deltas and deriva- 
tives of deltas) which appear when the long distance po- 
tential U\,{r) is computed in perturbation theory. These 
terms are in fact ambiguous (and hence unphysical) and 
depend on the regularization scheme used in the pcrtur- 
bative calculation, but are not essential since they do not 
contribute to the absortive part and hence to the cor- 
responding spectral function 0. The important point 
is that these ambiguous distributional terms never con- 
tribute to the long range part provided r c > since in 
this case a compact support for distributional terms is 
guaranteed. This is a clear advantage of the radial cut- 
off we are using. In fact, our main motivation for carrying 
out the analysis in coordinate space is this clean separa- 
tion between short and long distances. Many regulators, 
mainly those in momentum space, do not fulfill this con- 
dition, since the regulator effectively smears these short 
distance terms and intrudes somewhat into the long dis- 
tance region 3 . 

Finally, we need a matching condition connecting the 
inner and outer regions, 

«fc, S ( r c) 

1 — T = i \ 6 

Uk,L(r c ) Uk,s\rc) 

Viewed from the outer region this relation corresponds to 
an energy dependent boundary condition at a given short 
distance cut-off radius, r c . Because of elastic unitarity we 
expect the state to be normalized, so that if we use a box 
of size a as an infrared regulator, we have 

J ° u kiS (r) 2 dr + J u kih (r) 2 dr = 1 (7) 

with a much larger than the range of the potential. 
This equation gives a quantitative separation between 
long distance and known physics and short distance and 
unknown physics. Obviously, any effective description 
based on the long range part should fulfill 

/ C u k , s (r) 2 dr -C I (8) 
Jo 



which corresponds to the requirement of a small wave 
function in the inner unknown region. This is our ba- 
sic condition from which most of our results follow. It 
should be realized that here long and short distances are 
intertwined through the matching condition, Eq. In 
particular, an arbitrarily growing function at the origin 
cannot fulfill this condition even if value of the short dis- 
tance cut-off is taken. 

This kind of pathological situation actually occurrs 
when dealing with the deuteron channel in the theory 
with no explicit pions with OPE potential j^l) and 
TPE potential at NLO in the Weinberg counting (see 
Sect. IVliTl and IVIII CI below - ) . As an instructive and en- 
lightening example, we illustrate the situation in Fig. ^ 
in the deuteron state for the quantity 

pro 

P{r c ) = / {u{rf + w{r) 2 )dr (9) 
Jo 

(for notation see Sect. El using the LO, NLO and NNLO 
potentials, Eq. (JTJ, in the outer region, r > r c , and 
matching to a free particle in the inner region r < r c 
(the precise form of the wave function inside turns out 
not to be essential) 4 . As one sees, there are cases such 
as a pure short distance theory where the asymptotic 
D/S ratio, rj, and the deuteron binding energy are fixed 
to their experimental value. For r c < 1.4fm one has a 
miminum probability of about 20% in the inner region, 
and then P(r c ) starts to grow. In this case the description 
cannot be considered effective below that critical value. 
A similar situation occurs in the NLO TPE case where 
one has a decreasing inside probability until one reaches 
r c ~ 0.8fm where it takes its minimum value, about 7% 
and then starts increasing. In contrast, for the LO (OPE) 
case where the deuteron binding energy is fixed to the ex- 
perimental one and 77 is predicted |41j to have the value 
77 = 0.02633 from the regularity condition of the wave 
function at the origin and also the NNLO (TPE) case 
where both numbers are fixed to experiment, the prob- 
ability in the interior region is controlled and generally 
decreasing. This is a first and transparent illustration 
that the description based on any preconceived power 
counting is not necessarily consistent with the fact that 
short distance ambiguities are under control, since the 
wave function in the inner and unknown region does not 
become arbitrarily small as the cut-off is removed. An- 
other possibility is to keep a finite cut-off, and we will 
analyze this point in more detail in Sect. IV111 CI 

As we will discuss in the next Sect. IIIII tight con- 
straints on the structure of short distance counterterms 
must be fulfilled if the requirement P(r c ) — ► as r c — > + 
is imposed. 



For in stance, the widely used gaussian regulator in momentum 
space |2a. |29L I42l of a local potential in cordinate space corre- 
sponds indeed to a convolution of the original potential smeared 
over the region of size 1/A. 



4 In a momentum space formulation this is somewhat equivalent 
to cut-off the Lippmann-Schwinger equation above a given value 
A. 
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FIG. 1: Probability inside the short distance region as a func- 
tion of the short distance cut-off radius r c for the deuteron 
state for different approximations to the potential. In all cases 
the deuteron binding energy is fixed to the experimental value. 
— tv means a pionless theory where in addition n = 0.0256 is 
also fixed to the experimental number. OPE means our pre- 
vious work |4l| where r\ = 0.02633 is deduced from the regu- 
larity condition of the wave function at the origin. OPE + r\ 
coresponds to the OPE case when n = 0.0256 (the experimen- 
tal value). TPE NLO and TPE NNLO (Set IV) correspond 
to the higher terms in Eq. Q where also the experimental n 
value is taken. 



Wigner bounds on the short distance 
contributions 



To proceed further, we derive with respect to the en- 
ergy Eq. JSJ) and one immediately gets after some algebra, 



d 
dk 2 



ru 'k,s( r c) 
Uk,s(r c ) 



Jo" u k ^{r) 2 dr 



< 



(10) 



whenever r c is not a zero of the wave function. 

The short range theory can be characterized by an ac- 
cumulated phase shift Ss(k, r c ), given by the solution of 
the truncated short range problem, 

«i,sW = sm(kr + 8 s (k,r c )) , r > r c (11) 

which fulfills, from Eq. l(T0")l. 

-L[k cot (kr c + S s (k,r c ))} <0 (12) 

a condition equivalent to Wigner's causality condi- 
tion [34| . Using an effective range expansion for the short 
distance phase shift 



k cot 5 s 



1 



1 



-Jo,sk 



we get the Wigner bound for the effective range 



ro,s < 2r c 



1 



r c 



(13) 



(14) 



where ao,s and ro.s represent the scattering length and 
effective range when the short distance potential Ug{r, r') 
is switched on from the origin up to the scale r c or equiv- 
alently when the long distance potential is switched off 
from infinity down to r c (see Refs. 39, 40] for more de- 
tails). In a theory where the long distant potential is 
absent Uh(r) — 0, i.e., a pure short distance descrip- 
tion, we have the obvious result that the short distance 
threshold parameters coincide with the physical param- 
eters ao.s = cto and r Q $ = Tq. Thus, 



r < 2r c 



1 - 



00 



' c 

3a 



, U L (r) = 



(15) 



which implies that r$ < for r c — ► 0. With the exper- 
imental values one gets the lowest short distance cut-off 
compatible with causality to be r c = 1.4fm. For a given 
long distance potential we just solve the equations from 
infinity inwards and look for the point where the Wigner 
condition is first violated. In Fig. [21 we plot the evolution 
of the Wigner bound on the effective range for the differ- 
ent approximations to the potential according to the ex- 
pansion (JIJ as a function of the short distance cut-off ra- 
dius. As we see, the lower bound on the radius is pushed 
towards the origin, and in fact for the NNLO approxima- 
tion there is no lower bound at all. Thus, only for the 
NNLO TPE potential can one build the full strength of 
the experimental effective range without violation of the 
Wigner condition. We will see more on this in Sect. HVl 
So far the discussion has been carried out for a fixed 
value of r c . If we change the short distance radius, r c — > 
r c + Ar c , we can use the matching condition, Eq. 10, 
to evaluate the change seen from the outer region. This 
results into a variable phase equation which has been 
analyzed extensively in our previous works |39l Eoj . If 
we take the limit r c — > 0, we get that the outer wave 
functions fulfills 



d 
dk 2 



<(0 + ) 
u k {0+) 



= 



(16) 



where the label L has been supressed. Thus, the bound- 
ary condition becomes energy independent when the limit 
r c — * + is taken if the inner wave function becomes ar- 
bitrarily small. Note that the limit is taken from above 
such that r c > 5 . As we have said already, this jus- 
tifies not considering contact terms in the potential. A 
direct consequence is that we get also that with the ex- 
ception of the short distance scattering length as,o(0 + ), 
which can be fixed by some renormalization condition 



To illustrate this point, let us note that for potentials which do 
not diverge too strongly at the origin r' 2 U(r) — » such as OPE 
in the 1 Sq channel, there is some irreversibility in the process of 
integrating from exactly r c = out (which requires the regular 
solution, u(0) = 0) or integrating in towards the origin from 
above r — » + (which generally invo lves the irregular solution, 
u(0) 7^ 0). See the discussion in Ref. I'jll and in Sect. ITTTT?! 
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FIG. 2: Wigner bound on the effective range of the short dis- 
tance interaction, (r.h.s. of Eq. 1141 ') as a function of the short 
distance cut-off radius for the contact interaction (ERE), 
the OPE interaction (LO), the TPE (NLO) and NNLO. We 
use the conditions ao — — 23.74fm and ro = 2.77fm when 
r c l/m n . Causality is violated when the curve crosses 
the experimental value. Straight horizontal lines correspond 
to the values obtained for the effective range when the short 
distance initial condition ro,s(0 + ) = is assumed. 



(see Sect. IHIfl . the remaining short distance threshold 
parameters are also zero in this limit, 



distance properties of chiral potentials have nothing to 
do with short distance properties of the "true" potential, 
but renormalization and finiteness requires a very pre- 
cise behaviour of the wave function when approaching 
the origin from long distances. In this section we classify 
the undetermined constants depending on the attractive 
or repulsive nature of the corresponding potentials in the 
single channel case and the eigenvalues of the potential 
matrix in the coupled channel case. In coordinate space 
and disregarding relativistic corrections the potentials in 
Eq. are local and energy independent 6 . An important 
condition on the short distance behavior of the wave func- 
tions are the orthogonality constraints between states of 
different energy. For a regular energy independent po- 
tential these constraints are automatically satisfied, but 
for singular potentials they generate new relations rel- 
evant to the NN interaction. Our approach is not the 
conventional one of adding short distance counterterms 
following an a priori power counting on the long dis- 
tance potential. Rather it is the power counting in the 
potential what uniquely determines the admissible form 
of the short distance physics if we want to reach a finite 
limit when the regulator is removed. This can only be 
achieved by choosing the regular solution at the origin, 
i.e. u(0) — 0. In Sect. IV 111 CI we will show that irreg- 
ular solutions generate divergent results after renormal- 
ization. Although this may look as a potential drawback 
of removing the cut-off it provides valuable insight on the 
form of the potential and the validity of the expansion 

(see sect. iynrni . 



ro,s(0 + ) = f 2 ,s(0 + ) = 0, .. 



(17) 



This energy independence of the boundary condition at 
the origin insures the orthogonality conditions between 
different energy states, 



o+ 



Uk(r)u k > (r)dr = S(k - k') . 



(18) 



A. One channel case 

Let us first review the single channel case in a way that 
results for the coupled channel situation can be easily 
stated. The reduced Schrodinger equation for angular 
momentum I is 



U(r)u + = k 2 i 



(20) 



III. SHORT DISTANCE BEHAVIOR OF 
CHIRAL POTENTIALS, ORTHOGONALITY 
CONSTRAINTS AND THE NUMBER OF 
INDEPENDENT CONSTANTS 

As we have said, chiral NN potentials, Eq. (J), al- 
though decay exponentially at large distances, become 
singular at short distances, where one has 



U(r) 



MC n . 2 
— — (1 + air + a 2 r 



(19) 



To avoid any misconception let us emphasize that the 
short distance behaviour of a long distance potential 
should be regarded as a long distance feature, i.e. a long 
wavelength property, since different long distance poten- 
tials yield different short distance behaviours. The short 



For a power law singular potential at the origin of the 
form U(r) = MC n /r n = ±{R/r) n /R 2 with n > 2 and R 



6 In momentum space and up to NNLO the long distance part of 
the potential depends on the momentum transfer q only and not 
on the total momentum k. Essential non-localities, i.e. contribu- 
tions of the form V(q, k) = L(q)k 2 with L(q) a non-polynomial 
function, depe nd weakly on the t ota l momentum and appear 
first at N 3 LO El B IS 13 E3| due to relativistic 1/M 2 
one loop contributions. In coordinate space this weak non- 
locality corresponds to a modification of the kinetic energy term 
in the form of a general self-adjoint Sturm-Liouville operator, 
— u"(r) — > — (p(r)w'(r))', with a singular p(r) function at the 
origin and exponentially decaying at long distances. The present 
formalism can in principle be extended to include these features, 
and will be discussed elsewhere. Nevertheless, according to the 
results of Sect. lVll fsee Table ITVL on the loss of predictive power 
already at NNLO there is a lack of phcnomcnological motivation. 
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the length scale dimens ion, th e de Broglie wavelength is 
given by l/k(r) = 1/ \/\U(r)\ and the applicability con- 
dition for the WKB approximation reads (l/fc(r))' <C 1, 
so that for distances r <C R(n/2) 2 ^ 2+n ^ one has a semi- 
classical wave function 0, 0, . Keeping the leading 
short distance behavior one gets for attractive and repul- 
sive singular potentials and any angular momentum the 
following behaviour for the regular solutions 



ux(r) - C A (L) 



n/4 



R 



n — 2 \ r 



for U A 
UR(r) 
for U R 



1 / R 
+ R? I 7 



»/4 



exp 



R 



n — 2 V r 



(21) 
(22) 

(23) 

(24) 



respectively. Here C 'a and Cr are normalization con- 
stants and ip an arbitrary short distance phase. In the 
repulsive case we have discarded the irregular solution 
(a similar exponential with a positive sign) which would 
not allow to normalize states. For an attractive singu- 
lar potential there is a short distance unknown parame- 
ter. This phase could, in principle, be energy dependent. 
Chiral potentials are, however, local and energy indepen- 
dent at NNLO at all distances, and become genuinely 
energy dependent at N 3 LO due to relativistic 1/M 2 cor- 
rections [3 EH H [H Ell 7 ■ Thus, if we require or- 
thogonality of states with different energy (positive or 
negative) we get 



= 



L k LL p 



j, p a k 



R 



sin (ip(k) - ip(p)) 



(25) 



Hence, the phase ip is energy independent and could be 
fixed by matching the solution to the asymptotic large 
distance region (we assume a short range potential), e.g., 
by requiring a given value of the scattering length, a;, at 
zero energy. In this way, a new and physical scale appears 
into the problem which is not specified by the potential. 
This is equivalent to the well known phenomenon of di- 
mensional transmutation. Another possibility is to fix ip 
from a given bound state energy, E = —B. The new scale 
entering the problem is the corresponding wave number, 
7 = \J MB. Note that although neither ai nor 7 can 



7 The subthreshold energy dependence from the virtual pion pro- 
duction channel NN — > NNtt which is in principle N 3 LO, disap- 
pears since in the heavy baryon limit the threshold s^jviv = 
(2M + m) 2 = 4(M 2 + k 2 ) translates into a CM momentum 
k = v / m(M + m/4) — > 00. 



be predicted from a singular potential, the orthogonal- 
ity constraint does predict a correlation between them 
through the potential. Likewise, the phase shifts 5i can 
be deduced from either ai or 7 by taking the same short 
distance phase ip. In the repulsive case there is no dimen- 
sional transmutation since the orthogonality condition 
follows from regularity at the origin, and the potential 
fully specifies the wave function. In this case, the scatter- 
ing length and the spectrum are completely determined 
from the potential as for standard regular potentials. 



B. Coupled channel case 

We turn now to the two coupled channel case where 
the wave functions are denoted by a column vector (u, w) 
(for some particular cases see e.g. Refs. HHEiia). If 
we assume that at short distances the reduced potential 
behaves as 



U->M 



C n 



(26) 



where C„ is a symmetric matrix of Van der Waals coef- 
ficients. Diagonalizing the matrix C n we get 



cos a 
- sin( 



sin f 

cos( 



C n ,+ 
Cn,. 



cos V — sin t 
sin 8 cos 8 



(27) 



where C n ,± are the corresponding eigenvalues and 8 the 
mixing angle. Thus, at short distances we can decouple 
the equations to get 



cos 8 sin 8 
- sin 8 cos 8 



(28) 



where (u+, ti_) are regular solutions as in the single chan- 
nel case. So, in the two channel situation we have three 
possible cases depending upon the sign of the eigenvalues. 



both eigenpo- 

-2 



1. Both eigenvalues are negative, i.e. 

tentials are attractive and MC n .+ = — i?™ - ^ and 
MC nj _ = — R™T 2 with R± the corresponding scale 
dimension. In this case the short distance eigen- 
solutions are oscillatory and there are two undeter- 
mined short distance phases, tp+ and tp-. Moreover 
for two states, (uk,Wk) and (u p ,Wp), with different 
energies we get the orthogonality constraint 



= u' k u p — u p u' k + w' k w p — w p w' k 



R. 



sin (ip+(k) - <p+{p)) 



1 

rZ 



sin (¥>-(fc) - ¥-(p)) ■ 
(29) 



2. One eigenvalue is negative and the other is posi- 
tive, MC„ !+ = R\- 2 and MC n _ = -R^ 2 . One 
short distance eigensolution is a decreasing expo- 
nential and the other is oscillatory, so we have one 
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short distance phase (p. In this case for two states 
(itfc, w k ) and (u p , w p ) with different energies we get 
the orthogonality constraint 

= u' k u p — u p u' k + w' k w p — w p w' k 

= ^-sm(<p + (k)-cp + (p)) . (30) 

3. Both eigenvalues are positive, MC„, + = i?"~ 2 and 
MC n ,~ = i?"~ 2 . Then, both short distance eigen- 
solutions are decreasing exponentials. There are no 
short distance phases. In this case the orthogonal- 
ity relations are automatically satisfied. 

This simple argument can be easily generalized to any 
number of coupled channels. The number of undeter- 
mined short distance phases corresponds to the number 
of attractive eigenpotentials at short distances. Orthog- 
onality of the wave functions requires that all these short 
distance phases fulfill a generalized condition of the form 
of Eq. {EJ. 

The orthogonality conditions require the determina- 
tion of the short distance phases, as we did in Ref. |4l| 
for the OPE case. This requires in general an improve- 
ment on the short distance behaviour to high orders. An 
alternative method is to impose the orthogonality con- 
straints either in the single or coupled channel case by 
integrating in from infinity for a fixed energy, either pos- 
itive or negative, and then impose the condition at a 
sufficiently short distance cut-off radius r = r c . In the 
single channel case one would get the condition, 

Uk{r c ) u Q (r c ) 

if the zero energy state is taken as the reference state. 
An analogous relation holds for the coupled channel sit- 
uation, namely 

= u k (r c )u a (r c )' - u k (r c )'u (r c ) 

+ w k (r c )wo(r c y - w k (r c )'w (r c ) . (32) 

Obviously, in this procedure cut-off independence must 
be checked. For the TPE chiral potentials analyzed in 
this paper we find that r c = 0.1 — 0.2fm proves a suffi- 
ciently small value of the short distance cut-off. 

C. Power Counting, counterterms and short 
distance parameters 

As we see, the number of independent parameters is 
determined from the potential, although their value can 
be fixed arbitrarily, by some renormalization condition 
like, e.g., fixing scattering lengths to their physical value. 
This removes the cut-off in a way that short distances 
become less and less important. Now, if the potential 
is regular, i.e., r 2 |L7(r)| — > 0, one may choose between 



the regular and irregular solution 8 . In the first case the 
scattering length is predicted while in the second case the 
scattering length becomes an input of the calculation. In 
either case the wave function is still normalizable at the 
origin. Singular potentials at the origin, i.e. fulfilling, 
r 2 |[/(r)| — ► oo, do not allow this choice if one insists on 
normalizability of the wave function at the origin. If the 
potential is repulsive, the scattering length is fixed while 
for an attractive potential the scattering length must be 
an input parameter. Furthermore, orthogonality of dif- 
ferent energy solutions requires an energy independence 
of the boundary condition, so that in all cases the effec- 
tive range, and higher order threshold parameters cannot 
be taken as independent parameters, in addition to the 
scattering lengths. 

This can be translated into the language of countert- 
erms quite straightforwardly. In momentum space, fix- 
ing ao arbitrarily corresponds to take a constant Co cut- 
off dependent and energy independent contribution to 
the potential Vo(k', k) in the Lippmann-Schwinger equa- 
tion. Likewise, fixing ro can be mapped as adding a term 
Cii}i 2 + k' 2 ) to the potential. For higher coupled chan- 
nel partial waves one fixes the scattering length otiji one 
has instead terms of the form Ci> ik a k l in the potential 
V Vt i(k',k). 

The OPE potential in the singlet 1 So is regular at 
the origin and hence one can take On as an indepen- 
dent parameter or not (see Refs. |3|J,|4(|.) Actually, the 
smallness of the scattering length for the regular solution, 
suggests using the irregular solution. In the Weinberg's 
counting of the potential, at NLO one has TPE contri- 
butions in the potential. At short distances they be- 
have as an attractive 1/r 5 potential (see Sect. IVIIIf) . and 
then ao must be an independent parameter. At NNLO 
one has, again, a singular attractive 1/r 6 potential (see 
Sect. HVft . and thus an adjustable scattering length. This 
looks quite natural because increasing the order in the 
potential has a meaning and we can always compare the 
effect in the phase shifts of having a higher order poten- 
tial with the same scattering length (See Sect. IVIIIfl . 
In this construction, if the next term in the expansion 
turned out to be more singular and repulsive the scatter- 
ing length would be fully predicted from the potential. 

The OPE potential in the triplet 3 Si - 3 D x coupled 
channel corresponds to case 2) and hence one has one free 
parameter in addition to the OPE potential parameters. 
One may choose this parameter to be the deuteron bind- 
ing energy (or alternatively the triplet S-wave scattering 
length) . Any other bound state or scattering observables 
are predicted. This case was treated in great detail in our 
previous work |4l| 9 . In the NLO TPE potential we have 



Both cases comply to the normalizability condition at the origin, 
Eq. JHJ 

Relevant previous work on this channel was also presented in 
Ref. l38l and |4(J where the orthogonality conditions where not 
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case 3) because both eigen potentials present a repulsive 
1/r 5 singularity (see Sect. I Vlllj) and one would predict 
all observables from the potential parameters. Finally, in 
the NNLO TPE potential we have case 1) corresponding 
to an attractive- attractive (see Sect. [Vj and two addi- 
tional parameters need to be specified for a state with a 
given energy. The orthogonality condition imposes a re- 
lation between two states of different energy, so that for 
all energies in the triplet channel one has three indepen- 
dent parameters. We will take these three parameters 
to be the deuteron binding energy, the asymptotic D/S 
ratio and the S— wave scattering length. The trend one 
observes when going from LO to NNLO is quite natural; 
as usual in ChPT one has more parameters at any order 
of the approximation. The NLO approximation poses, 
however, a problem since one seems to have more predic- 
tive power than at LO (See Sect. IVf 1 1 1 for more details 
on the consequences of using our renormalization ideas 
literally for the conventional NLO potential). 

For the NNLO TPE triplet 3 Si - 3 D 1 channel we fix 
the deuteron binding energy, or equivalently 7, and the 
asymptotic D/S ratio rj by their experimental values. 
This fixes the short distance phases ^+(7) and (p-("/). 
Next, if we use an a or [3 (see below for a definition) zero 
energy scattering state we have in principle two short dis- 
tance phases (<^ a , + (0), <^q, ; _(0)) and (<^+(0), <^/3,_(0)) 
which can be related to the (ao,o;o2) and (ao2,o;2) scat- 
tering lengths respectively. Using the orthogonality con- 
straints to the deuteron bound state one can then elimi- 
nate CH02 and 012 and treat ao as a free parameter. Thus, 
in the triplet 3 Si — 3 D\ channel we can treat 7, r] and 
aQ as independent parameters. Once these parameters 
have been fixed we can actually predict the correspond- 
ing phase shifts since any positive energy state must be 
orthogonal both to the deuteron bound state and the 
zero energy scattering states. This result is a direct con- 
sequence of the singular Van der Waals attractive behav- 
ior of the TPE potential at the origin. It is remarkable 
that this same set of independent parameters was also 
adopted in Ref. |54| within the realistic potential model 
treatment. 

Conflicts between naive dimensional power counting 
and renormalization have been reported recently already 
at the LO (OPE) level 01 where it is shown that even 
the 3 Po partial wave depends strongly on the cut-off in 
momentum space (a gaussian regulator is used) if accord- 
ing to the standard Weinberg counting no counterterm is 
added. The requirement of renormalizability makes the 
promotion of one counter term unavoidable for channels 
which present an attractive singularity. This promotion 
is the minimal possible one compatible with finiteness, 
because in a coupled channel problem one could think 
in general of three counterterms. From this viewpoint 
the choice of just one counterterm in triplet channels is a 



bit mysterious. Our discussion in coordinate space agrees 
with these authors in the OPE potential, and actually al- 
lows to identify a priori the necessarily promotable coun- 
terterms as non trivial boundary conditions at the ori- 
gin for singular attractive potentials. Moreover, we also 
see that the promotion of only one counterterm in the 
triplet channels with an attractive-repulsive singularity 
invoked in Ref. |42j is also maximal, since any additional 
counterterm would also produce divergent results (see 
Sect. lVIII Cj) . Thus, we see that although power counting 
determines the long distance potential, the short distance 
singular character of the potential does not allow to fix 
the counterterms arbitrarily. 

To conclude this discussion, let us mention that the 
short distance phases, whenever they become relevant 
play the role of some dimensionless constants which de- 
pend exclusively on the form of the potential, but not on 
the potential parameters . For the same reason they 
can be taken to be zeroth order in the power counting 
used to generate the chiral potential in Eq. although 
they are subjected in general to higher order corrections. 
In this sense, the form of the short distance interaction 
is dictated by the potential only, and cannot be consid- 
ered as independent information (See also Ref. [5^ for a 
similar view on the three body problem in the absence of 
long distance potentials). 



IV. THE SINGLET 'So-CHANNEL 

For the singlet 1 S'q channel one has to solve 

-u"{r) + Ui Sa (r)u(r) = k 2 u{r) (33) 

At short distances the NNLO NN chiral potential behaves 
as HlHd 



3 5 2 



128/ 4 7r 2 r 6 
_R*_ 

7-6 



(-4 + 15.g 2 + 24c 3 - 8c 4 ) 



(34) 



which is a Van der Waals type interaction with typical 
length scale R = (-MC e ) 1/4 . Here, c, = Mc l . The value 
of the coefficient is negative for the four parameter sets of 
Table[U so the solution at short distances is of oscillatory 
type, Eq. <|2*T|l with n — 6, and 



u(r) 



1 /R 

2 17 



■If 



(35) 



considered. See the discussion at the end of Sect. V in Ref. |4 ll 



where there is a undetermined energy independent phase, 
ip, and A is a normalization constant. Note that the 
corresponding Van der Waals radius is quite sensitive to 
the choice of chiral parameters. 
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TABLE I: Short distance Van der Waals coefficients for the NNLO chirai potential in singlet 1 So and the 3 Si — 3 D\ triplet 
channels. 



Set 


Source 


c^GeV" 1 ) 


c 3 (GeV- 1 ) 


c 4 (GeV- 1 ) 


MC 6 (fm 4 ) 


MC 6 ,+ (fm 4 ) 


MC 6 ,- (fm 4 ) 


9 (degrees) 


Set I 


(BM) nN [31] 


-0.81± 0.15 


-4.69± 1.34 


3.40 ± 0.04 


-8.74 


-16.96 


-5.63 


140.8 


Set II 


(RTdS) NN [11] 


-0.76 


-5.08 


4.70 


-10.19 


-21.45 


-5.58 


170.0 


Set III 


(EMa) NN [20] 


-0.81 


-3.40 


3.40 


-6.45 


-14.68 


-3.35 


140.7 


Set IV 


(EMb) NN [28] 


-0.81 


-3.20 


5.40 


-7.28 


-20.18 


-1.86 


182.6 



A. Low energy parameters 

For the zero energy state we use the asymptotic nor- 
malization at large distances 



a 



Then, the effective range is given by 



'o 



2 / dr 
'o 



on 



u (rf 



(36) 



(37) 



We can use the superposition principle for boundary con- 
ditions 



u (r) = uo,c(r) - —u , s (r) 
ao 



(38) 



where uq }C (t) — > 1 and Uo, s (t) —* r correspond to cases 
where the scattering length is either infinity or zero re- 
spectively. Using this decomposition one gets 



ro 



= A 



B 

a 



C 



where 



.4 
D 

C 



M 1 ~ u l.c) 



-4 / dr(r - u , c u 0iS ) , 
Jo 

/>oo 

2 / dr(r 2 - «g ) , 
Jo 



(39) 

(40) 
(41) 
(42) 



depend on the potential parameters only. The interest- 
ing thing is that all explicit dependence on the scattering 
length a is displayed by Eq. J2H ) . In a sense this is the 
non-perturbative universal form of a low energy theorem, 
which applies to any potential regular or singular at the 
origin which decays faster than a certain power of r at 
large distances (for an analytical example with the pure 
Van der Waals potential U = -R 4 /r 6 see Sect. NTfy . 
Since the potential is known accurately at long distances 
we can visualize Eq. (|39|l as a long distance correlation 
between ro and ao. Naturally, if there is scale separa- 
tion between the different contributions in the potential, 
Eq. (JTJ, we expect the coefficients A,B and C to dis- 
play a converging pattern. This is exactly what happens 



(see Eq. (|43|) and Eq. (|89|l below) although not compat- 
ible with a naive and perturbative power counting (see 
Sect.llXl. 

In the 1 S'o the TPE potential becomes singular and at- 
tractive at short distances. Nevertheless, already at this 
point one can see the dramatic difference between attrac- 
tive and repulsive singular potentials. In the attractive 
case the short distance phase allows to choose the scat- 
tering length independently on the potential, hence the 
coefficients A,B and C are uncorrelated with ao- For a 
singular repulsive potential, however, ao as well as A,B 
and C are determined by the potential. If one assumes 
A,B and C to be independent on ao in the repulsive 
case, this can only be possible due to an admixture of 
both the regular and irregular solutions, the latter will 
dominate at short distances and the effective range will 
diverge ro — > — oo. This fact will become relevant in 

Sect.|vnfni 

Obviously, for the chirai TPE potential, Eq. (QJ, the 
coefficients have to be evaluated by numerical means and 
they are finite. We expect that these coefficients scale 
with the relevant scale of the potential. If long distances 
dominate A ~ 1/m. B ~ 1/m 2 and C ~ 1/m 3 but then 
ro ~ 1/m. On the contrary if short distances dominate 
A ~ R, B - R 2 and C ~ R 3 and r Q - R. The real situa- 
tion is somewhat in between, but it is clear that A is far 
more sensitive to short distances than C. Actually, for a 
large scattering length, as it is the case in the 1 5o channel, 
the coefficient A dominates. Note that unlike the stan- 
dard approaches, where a short distance contribution to 
the effective range is allowed (in the form of a momentum 
dependent counterterm Ci{k 2 + k' 2 ) ), we build ro solely 
from the potential and the scattering length ao- This is a 
direct consequence of the orthogonality relations, which 
preclude energy dependent boundary conditions for the 
local and energy independent chirai TPE potential. 

In Fig. (J3J we show the dependence of the effective 
range as a function of the short distance cut-off radius r c , 
i.e. replacing the lower limit of integration in Eq. I|37|) . for 
values between 2 and O.lfm and taking the experimental 
value of the scattering length ao = — 23.74fm. As we see, 
the short distance behaviour is well under control and 
nicely convergent towards the experimental value. This 
dependence also illustrates that an error estimate based 
on varying the cut-off between certain range is only a 
measure on the size of finite cut-off effects, rather than 
a measure on the error. The linear behaviour observed 
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TABLE II: Threshold parameters in the singlet 1 So chan- 
nel for the different sets of parameters ci, C3 and C4 given 
in Table |T] We compare our renormalized results given by 
the cut-off independent universal formula 13911 for ro and its 
extension for V2 to finite cut-off NN calculations using their 
scattering length as an input. The difference is attributed to 
finite cut-off effects. 
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-0.23 


Set III 
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Set IV 


This work 


oo 


Input 
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Nijm II 


[43., 44] 
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2.67 


-0.48 
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[43,44] 
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-0.49 
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-23.74(2) 
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FIG. 3: Effective range for the TPE potential in the singlet 
1 So channel when the lower limit of integration in Eq. 1371 
is taken to be a short distance cut-off and the experimental 
value of the scattering length ao = — 23.74fm is taken. We 
show Sets I, II,III and IV (see main text). 



at small r c is a consequence of the dominance of the first 
term in Eq. (|37[) as compared to the second term where 
the wave function contribution vanishes as ~ r\. Let us 
remind that in the conventional treatments a countert- 
erm C% is added to provide a short distance contribution 
to the effective range parameter and the result is fitted 
to experiment so that ro becomes an input of the cal- 
culation. Obviously one expects Ci to depend on the 
rcgularization scale. As shown in Fig. © the size of the 
counterterm C% at r c — > must be numerically small, 
since the TPE potential provides the bulk of the con- 
tribution. This agrees with the orthogonality constraint 
which requires C2 — > when r c — > 0. 

Numerically we find the following renormalized rela- 
tions in the singlet channel for the OPE and NNLO TPE, 



ro = 1.308062- 



ro 



2.670963 



r = 2.715075- 



2.586862 



2.616830 



4.547741 



5.192606 


a 


2 

"o 


5.755234 


6.031119 




ao 


2 


5.847358 


6.093430 




ao 


5 


5.584383 



5.916900 


a Q 


2 

"o 


5.640921 


5.952694 




a Q 


2 

"0 



(OPE) 
(Setl), 



(Set II), 

(Set III) , 

(Set IV) 

(43) 

As we see, the coefficient dependent of cvq 2 is not very 
sensitive to the choice of the coefficients c l5 C3, C4 and the 
OPE potential already provides the bulk of the contribu- 
tion. On the other hand, the coefficient independent on 



cvo changes dramatically when going from OPE to TPE, 
suggesting that the effect is clearly non-perturbative. A 
direct inspection of the integrands for the A,B and C co- 
efficients shows that A picks its main contribution from 
the short distance region around 1 fm, whereas for B and 
C the most important contribution is located around 3 
fm. One expects that different choices of coefficients C3 
and C4 influence mostly the A coefficient. Wc confirm 
this expectation analytically by only keeping the Van der 
Waals contribution to the full potential in Sect. IV 111 Wc 
emphasize, again, that A, B and C are intrinsic informa- 
tion of the potential; these values for the effective range 
stem solely from the NNLO chiral potential and the scat- 
tering length ao, without any additional short distance 
contribution. The closeness of these numbers to the ex- 
perimental value suggests that there is perhaps no need 
to make the boundary condition energy dependent if the 
cut-off is indeed removed, and that the missing O.lfm 
contribution can be clearly attributed to N 3 LO contri- 
butions in the potential. 

The results are summarized in Table [H] For pn 
we have the experimental values ao = —23.74(2) and 
r = 2.77(5). The previous formula, Eq. i|3"$j) yields 
r = 2.92,2.97,2.83,2.87 for Sets 1,11, III and IV re- 
spectively, which show a systematic discrepancy with the 
published values in several works (see References at the 
Table HT)l and also a systematic trend to discrepancy with 
respect to the experimental value. Our renormalized val- 
ues are always larger than the finite cut-off results. This 
seems natural since finite cut-off corrections diminish the 
integration region. Note also that the size of the discrep- 
ancy is larger than the experimental uncertainties and 
hence is statistically significant. 

The value of the effective range was not given in the 
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coordinate space calculation of Ref. 0] but the quality 
of the fit suggests that they get a value very close to the 
experimental one, ro = 2.75. The contribution to the 
effective range from the origin to O.lfm is about 0.2. In 
Ref. ^3 the cut-off is in coordinate space and an energy 
dependent boundary condition is considered. This means 
in practice cutting-off the lower integration in Eq. I|37|l 
at a = 1.4fm and adding a short distance contribution rs 
as to reproduce the experimental value. This introduces 
a new potential independent parameter. As we have ar- 
gued, in the limit a — > 0, the short distance contribution 
of the effective range should go to zero, as implied by the 
orthogonality constraints. For finite a, the orthogonal- 
ity constraint does not imply a vanishing short distance 
contribution to the effective range. 

For Set IV one could reach the upper experimental 
value by flipping the sign of c\ and keeping C3 and C4 
unchanged. For c\ = 2.43GeV~ 1 one gets ro = 2.78fm. 
The full experimental range would be covered by letting 
0.81GeV _1 < ci < 4.90GeV~ 1 . This is in total contra- 
diction to the expectations of nN scattering studies [3lJ . 
The insensitivity of our results with respect to the c\ co- 
efficient has to do with the fact that c\ only enters in 
the potential at short distances at order 1/r 4 which is 
sub-leading as compared to the leading Van der Waals 
singularity. This is another confirmation on the short 
distance dominance in the effective range parameter ro. 

Thus, according to our analysis, for the accepted values 
of chiral constants of Sets 1,11, III and IV used in previ- 
ous works, the difference in the value of the 1 S$ effective 
range could only be attributed to three pion exchange, 
relativistic effects and electromagnetic corrections. An- 
other possibility, of course, is to refit the chiral constants 
to our renormalized, cut-off free results. This will be 
discussed in Sect. IVII 



B. Phase Shift 

For a finite energy scattering state we solve for the 
chiral TPE potential with the normalization 

sin(/cr + 5q) 



u k (r) 



sin 5q 



(44) 



Again, if we use the superposition principle 

Uk(r) = u k , c (r) + k cot S u ktS (r) (45) 

with u k>c — > cos(fcr) and Uk, s ~ * sin(fcr)/fc and impose 
the orthogonality constraint with the zero energy state 
to get 

u 'k,c( a ^ + kcotS ° u 'kA a ) _ - a o u 'a,c( a ) + u 'oA a ) / 4g ) 
Uk,c(a) + k cot SoUk, s (a) -aou , c (a) +uo >s (a) 

Note that the dependence of the phase-shift on the scat- 
tering length is explicit; cot So is a bilinear rational map- 
ping of ao ■ Taking the limit a — > we get 

a A(k) - B{k) 




-20 
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FIG. 4: Renormalized phase shifts for the OPE and TPE po- 
tentials as a function of the CM np momentum k in the singlet 
1 S'o channel compared to the Nijmegen results |43j| for differ- 
ent parameter sets. The regular scattering wave functions for 
finite k are orthogonal to the zero energy wave functions. For 
the TPE potential we have taken the chiral couplings of set 
IV. 



whereas the functions A, B, C and T> are even functions 
of k which depend only on the potential and are given by 

A(k) = lim (u , c {a)u' k (a) - u' (a)u k , c (a)) 

a— *0 

B(k) = \hn(u k:C (a)u' 0s (a)~uo :S (a)u' kc (a)) 
C(k) = lim (u' (a)u k ,s(a) - u , c (a)u' k (a)) 

a— >0 ' ' 

V(k) = lim (uo iS (a)u' k (a) -u' (a)u kiS (a)) 

a—>0 



(48) 



The obvious conditions ^4(0) = D(0) = and B(0) = 
C(0) = 1 are satisfied. Expanding the expression for 
small k one gets the well known effective range expansion 



kcotS = -— + -r fc 2 + v 2 k 2 + 
ao 2 



(49) 



where v k is a polynomial in 1/ao of degree k + 1 . 

The renormalized phase shift is presented in Fig. 0] for 
Set IV. As we see the trend in the effective range ro and 
the t>2 parameter is reflected in the behavior of the phase 
shift 10 



k cot So 



a C(k)-V(k) 



(47) 



Let us mention that the momentum space calculation of Ref. |42| 
does not reproduce the physical and well measured scattering 
length. We have checked by an explicit solution of the Lipmann- 
Schwinger equation that the problem is more related to an in- 
sufficient number of Gauss points rather than to the value of the 
cut-off. 
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V. THE TRIPLET 3 Si - 3 D 1 CHANNEL 

The coupled channel 3 Si — 3 D\ set of equations read 
- u"(r) + Us Sl (r)u(r) + U El {r)w{r) = k 2 u{r) 



-w"(r) + U El (r)u{r) 



w(r) — k 2 w(r) , 
(50) 



At short distances the NN chiral NNLO potential behaves 
as UMM 



MC fi 3 



6, 3 Si 



MC, 



6 

6, -Ei 



r' 

MC, 



6 

6, 3 -Di 



(51) 



which is a coupled channels Van der Waals type interac- 
tion where the coefficients are given by 



MC aSl 

MC El 
MCs Dl 



3ff 2 



128/ 4 7T 2 

3V2g 2 



128/ 4 7T 2 

9ff 2 



(4 - Zg 2 + 24c 3 - 8c 4 ) 
(-4 + 3g 2 - 16c 4 ) 



(-l + 2.^ + 2c 3 -2c 4 ) 



32/ 4 7T 2 

If we diagonalize the corresponding matrix we get 



C6, 3 Si Ce.^i 



(52) 



cos 


sin 


— sin( 


? cos 


cos 6* 


— sin 


sin 


cos 6 




where C"6,± are the corresponding eigenvalues and 9 the 
mixing angle. They are listed in Table [I] for different pa- 
rameters choices of the chiral couplings c\ , c 3 and c 4 . We 
see that in all cases both eigenpotentials are attractive 
at short distances and hence the short distance behavior 
of the wave functions is of oscillatory type with n = 6. 
Defining the Van der Waals scales 



R± = (-MC 6 ,±) 1/4 
the short distance solutions read 



(54) 



cos 

— sin f cos 




Thus, we have two arbitrary short distance phases tp± for 
a given fixed energy which cannot be deduced from the 



u (T) = 0.0256) 
w (ti = 0.0256) 
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FIG. 5: The TPE deuteron wave functions, u and w, as a 
function of the distance (in fm ) compared to the Nijmegen II 
wave functions fill . . The asymptotic normalization u — ■» 
e ,r has been adopted and the asymptotic D/S ratio is taken 
r\ — 0.0256(4). We use the set IV of chiral couplings. 



potential and hence have to be treated as independent 
parameters. We will fix them to some physical observ- 
ables by integrating Eqs. (|50|l from infinity down to the 
origin. 



A. The deuteron 

In the deuteron k 2 = — -f 2 and we solve Eq. I|5(J|) to- 
gether with the asymptotic condition at infinity 



u(r) 
w(r) 



A 



se 



A D e~^ r 1 + 



jr ('yr)' 



(55) 



where 7 = V MB is the deuteron wave number, As is 
the normalization factor and the asymptotic D/S ratio 
parameter is defined by 77 = Ad /As- In what follows 
we use 7 and r\ as input parameters thus fixing the short 
distance phases ip± automatically. 

In this paper we compute the matter radius, which 
reads, 



(r 2 ) A f°° 

Y = l/ r 2 (u(r) 2 +w(r) 2 )dr 



(56) 



the quadrupole moment (without meson exchange cur- 
rents) 



1 f 00 

Qd= 20j r2u, M(2V2u(r) - w{r))dr . 



the deuteron inverse radius 



dr 



u(r) 2 + w(rY 



(57) 



(58) 
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FIG. 6: The dependence of the S-wave normalization As (in fm -1 / 2 , upper left panel) , the matter radius r m (in fm,upper 
right panel), the quadrupole moment Qd (in fm 2 , lower left panel), the D-state probability (lower mid panel) and the inverse 
radius (r" 1 } (in fm" 1 lower right panel) for the TPE potential on the short distance cut-off, r c (in fm) for Sets I,II,III and IV 
of low energy chiral constants. 



which appears in low energy pion-deuteron scattering, 
and the D-state probability 

POO 

Pd= w{rfdr. (59) 
Jo 

Following Ref. [41j we use the superposition principle of 
boundary conditions and write 

u(r) = us(r) + V u D(r) , 

w(r) = ws(r) + riwoir) , (60) 

where (us,ws) and (ud,wd) correspond to the bound- 
ary conditions at infinity, Eq. I|55|) . with As — 1 and 
Ad = and with As = and Ad — 1 respectively. Ob- 
viously, us,ud, ws and wd depend on the potential and 
the deuteron binding energy only, so that the dependence 
on the asymptotic D/S ratio rj can dc determined analyt- 
ically. The value is taken as a free parameter. The result- 
ing deuteron wave functions for Set IV are displayed in 
Fig. O and compared to the Nijmegen II results 0,01- 
One clearly sees the incommensurable ever increasing os- 
cillations already below r — 0.6fm. 

The short distance cut-off dependence of these 
deuteron properties using the experimental values for the 



deuteron binding energies and the asymptotic D/S ratio, 
r\ = 0.0256, can be looked up in Fig.HJ As one sees the 
cut-off dependence is well under control, so the infinite 
cut-off limit can be extracted without difficulty. 

Using the superposition principle of boundary condi- 
tions, Eq. I|60l ) the asymptotic S-wave ratio depends 
quadratically on 77 as follows 



I />oo />oo 

— = dr(u 2 s + w|) + 2rj dr{usu D + w s w D ) 
s Jo Jo 

/>OC 

+ rj 2 / dr(u 2 D +w D ). (61) 
Jo 



The coefficients of this second order polynomial depends 
on the potential and the deuteron binding energy. Sim- 
ilar relations hold for other observables. Evaluating the 
integrals numerically we get the following analytic corre- 
lations 
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Set I 

l/A% = 3.78888- 214.675 77 + 4489.43 ?7 2 

r 2 m /A 2 s = 5.47297- 54.1956 77 + 1295.89 rj 2 

Q d /A% = -0.342883 + 36.6449 77- 372.841 r? 2 

P D /A% = 2.10904- 184.824 77 + 4124.37 7i 2 

(r^ 1 )/^! = 3.58173 - 252.20 77 + 5186.53 t; 2 (62) 

Set II 

l/A 2 , = 3.01271 - 155.591 77 + 3363.94 t? 2 

rfjAl = 5.34737- 44.8896t;+ 1122.59 ?7 2 

Qd/A| = -0.296852 + 33.2406 77- 309.624 ?7 2 

Pd/A% = 1.44293 - 132.314 77+ 3098.89 t; 2 

(r^/Al = 2.52815- 171.36 77 + 3635.18 ?7 2 (63) 

Set III 

1/A 2 S = 4.65049 - 283.545 77 + 5902.53 r/ 2 

r 2 m /A 2 s = 5.58929 - 63.1854 77+ 1481.50 7i 2 

Qd/A% = -0.377779 + 39.2691 77- 420.250 ?7 2 

Pd/A% = 2.77521 - 241.491 77+ 5330.67 r/ 2 

(r'^/Al = 4.87639 - 356.21 77 + 7311.77 ?7 2 (64) 

Set IV 

l/A% = 3.40962- 190.713 77 + 4198.86 ?7 2 

r 2 n /A 2 s = 5.40066 - 49.2912 77+ 1232.81 ?7 2 

Q d /A% = -0.306469 + 33.9354 77- 318.598 t; 2 

P D /A% = 1.66525- 155.233 77 + 3681.89 ?7 2 

{r- l )/A 2 s = 3.14658- 226.20 77 + 4907.35 rj 2 (65) 

The numerical coefficients in these expressions depend 
on the deuteron binding energy and the TPE potential 
parameters, g, m, f, c\, C3 and a. The results for the 
deuteron properties are given in Table IIIII The uncer- 
tainties are due to changing the input 7 and 77 within 
their experimental uncertainties. We have checked that 
the short distance cut-offs a ~ 0.1 — 0.2fm generates much 
smaller uncertainties. The explicit dependence on 77 is 
displayed in Fig. [7] Again, we find a discrepancy in the 
case of Set I with the values quoted in the finite cut-off 
calculation. Remarkably, our renormalized results in co- 
ordinate space agree most with the momentum space cal- 
culation of Ref. 19] corresponding to Set IV. It is notice- 
able that this can be done without explicit knowledge of 
the counterterms used in that work in momentum space. 
This is precisely one of the points of renormalization; re- 
sults can be reproduced by just providing physical input 
data, and no particular reference to the method of solu- 
tion. Let us remind that the c\ , C3 and C4 were fixed from 
the perturbative study of NN peripheral waves where the 
cut-off sensitivity is rather small. Nevertheless, some sig- 
nificant discrepancies do also occur. 

For the parameter Set IV [2j| obtained by a N 3 LO 
fit to NN scattering data, our NNLO calculation almost 



reproduces exactly the numbers provided in that work. 
Furthermore, they turn out to be compatible with the 
experimental numbers at the la level within the uncer- 
tainty induced by the asymptotic D/S ratio 

One immediate lesson we learn from inspection of Ta- 
ble IIIII is that, regardless of the parameter set, only 
the experimental uncertainty in the asymptotic D/S ra- 
tio for the deuteron generates theoretical uncertainties 
about an order of magnitude larger then the experimen- 
tal ones. On top of this, one has also to take into ac- 
count other uncertainties, such as the one in g^NN and, 
of course, those induced by c±, C3 and C4, which gen- 
erally will generate larger uncertainties if all these pa- 
rameters are regarded as independent (see Sect. I VII be- 
low). In addition, there are systematic errors related to 
the accuracy of the expansion in the potential, Eq. {Tjl. 
In common with non-pert urbative finite cut-off calcula- 
tions dsmEJlillilln they are difficult t0 estimate 

a priori given the non-perturbative nature of our calcu- 
lation, but are bound to increase the error (see, however, 
our discussion in Sect. 11X1 below on non- integer power 
counting) . Given the insensitivity of our results with the 
short distance cut-off, the procedure used in Ref. [22I |23| 
of varying the cut-off becomes unsuitable in our case. 

For the deuteron channel one may conclude that the 
predictive power of the chiral expansion has reached a 
limit at NNLO. So, at present, we do not expect to make 
theoretical predictions in the deuteron to be more accu- 
rate than experiment. The inclusion of N 3 LO and higher 
orders may provide better central values but is unlikely 
to improve the situation regarding error estimates since 
new unknown coefficients in the potential appear and the 
induced uncertainties will generally increase. 

On the other hand, the slope for As and r m , Fig. 
suggests that it would be better to take the asymptotic 
S-wave normalization or the matter radius as input, since 
generated errors may be comparable or even smaller. For 
instance, if the matter radius r m is taken as input we get 
instead 77 = 0.0253(4) a compatible value with similar 
errors. However, if we take As — 0.8846(9) as input for 
Set IV we get 77 = 0.0255(1) a compatible value with 
the experimental one but with much smaller errors. The 
reduction of errors is also confirmed in Sets I, II and III, 
although the central values are a bit off. This result opens 
up the possibility of making a benchmark determination 
of the asymptotic D/S deuteron ratio from the chiral 
effective theory. Obviously, to do so, the chiral constants 
should be known with rather high accuracy, an illusory 



One may object that one should not use N 3 LO parameters to 
do a NNLO calculation, since they are obtained by fitting the 
same database. However, if there are finite cut-off effects the 
situation is not as clear. Finite cut-off eff ects are minimized 
in a N 3 LO calculation as shown in Ref. l2'l where the induced 
uncertainties are drastically reduced when going from NNLO to 
N 3 LO. Note that in our calculation there are no sizeable cut-off 
induced uncertainties already at NNLO. 
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expectation at the present moment. In this regard it 
would perhaps be profitable to pin down the errors for 
the chiral constants from peripheral waves. This point 
will be analyzed elsewhere 

Both the loss of predictive power and the very rare pos- 
sibility of making model independent theoretical predic- 
tions for purely hadronic processes using Chiral Pertur- 
bation Theory more accurate than experiment we seem 
to observe in low energy NN scattering is not new and 
has alre ady been documented for low energy ttit scatter- 
ing [43, ISO, |5l| and provides a further motivation to use 
chiral effective approaches. 



B. Low energy parameters 

The zero energy wave functions are taken asymptoti- 
cally as 12 



r 

1 , 

a 

r 



02_ oq2 \ _3_ _ r 3 
oto2 ct J r 2 15a 02 



(66) 



Using these zero energy solutions one can determine the 
effective range. The 3 Si effective range parameter is 
given by 



ro 



1- — I - u a (rf - w a (r) 2 
a 



dr . 



(67) 



In the zero energy case, the vanishing of the diverging 
exponentials at the origin imposes a condition on the a 
and P states which generate a correlation between ao , 
«02 and a.2. Using the superposition principle of bound- 
ary conditions we may write the solutions in such a way 

12 We correct an error in Eq. (45) of our previous work 14 ll where 
a.2 appears. The corrected numerical value is 02 = 6.693fm 5 . 
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(68) 



where the functions 1*1,2,3,4 and u>i,2,3,4 are indepen- 
dent on ao, ao 2 and a2 and fulfill suitable boundary con- 
ditions. The orthogonality constraints for the a and 
states read in this case 



u l u 0,f3 ~ U 7 U 0,/3 + Wy w 0,f3 ~ W 7 M 0,/3 



= 
= 



(69) 



yielding two relations between 7, ao2, ct2, tj and ao, 
meaning that two of them are not independent. Using 
the superposition principle decomposition of the bound 
state, Eq. llfiOfl . and for the zero energy states, Eq. 168(1 . 
we make the orthogonality relation explicit in ao, ao2, 
a2 and 77. If we would use ao, ao 2 , a2 as input parame- 
ters the orthogonality constraint is actually a non-linear 
eigenvalue problem for 7 and 77. The values of ao2 and 
a2 are not so well known although they have been deter- 
mined in potential models in our previous work |^. In 
contrast, 7, 77 and ao are well determined experimentally. 
Thus, in the deuteron scattering channel we will use 7, 
77 and ao as independent input parameters and ao2, oi 2 
as predictions. This same set of independent parameters 
was also adopted in Ref. | 54| within the high quality po- 
tential model treatment, although the role of the short 
distance Van der Waals singularity was not recognized. 
Fixing the experimental value of 7 we get the following 
relations for different parameter choices of c\, C3 and C4, 



Set I 

002 



Set II 

ao2 



2.01763 - 0.456461 a - 44.8947 r\ + 11.9351 a 77 
-0.314426+ 13.1555 7? 



-0.023522 + 1.04677 77- 11.6459 ?y 2 + a (0.008423 - 0.537856 r) + 9.39376 r/ 2 ) a 2 2 

a.2 = k 1 (70) 

a (-0.023901 + -qf a 



1.71745 - 0.373228 a - 33.4616 77 + 8.76639 a 77 
-0.228075 + 9.865911 77 



-0.030303 + 1.18083 77- 11.5032 ?7 2 + a (0.009559 - 0.566611 77 + 9.71850 7y 2 ) a% , . 

a 2 = i ^ — + — (71) 

a (-0.023118 + t?) 2 a 
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FIG. 7: The dependence of the S-wave normalization As (in fm -1 / 2 , upper left panel) , the matter radius r m (in fm,upper 
right panel), the quadrupole moment Qd (in fm 2 , lower left panel) (lower mid panel) and the inverse radius (r _1 ) (in fm" 1 
lower right panel) for the TPE potential as a function of the asymptotic D/S ratio r\. The boxes represent the experimental 
values. The predicted Quadrupole moments as well as the matter radius should be corrected for MEC's (accounting for adding 
O.Olfm 2 and 0.003fm respectively) on top of the potential result. We display the four sets of chiral coupling constants. 



Set III 

2.36659 - 0.550871 a - 59.1666 77 + 15.7488 a ry 



O-02 



-0.414962+ 17.2806 77 



-0.018755 + 0.937802 n- 11. 7229 rj 2 + a Q (0.007437 - 0.505434 r, + 9. 09925 rf) a L , s 

a<2 = — n ' '-L + -21 72 

a (-0.024013 + -qf a 

Set IV 

1.89526 - 0.418953a - 41.8369?/ + 10.8526007? 



ao2 



-0.279236+ 12.2978?7 



-0.023751 + 1.04857 Ty- 11.5733 ?7 2 + ao (0.008050 - 0.518604 7/ + 9.31798 t? 2 ) a \ 2 

Ot2 = 5 1 (73) 

a (-0.023118 + riY a 



The numerical coefficients appearing in these equations 
depend on the deuteron wave number 7 and the TPE 
parameters, g,f,m and ci, C3 and C4. The dependence 
on 77 for fixed values of ao within its experimental uncer- 
tainty is depicted in Fig. [SJ We see that for fixed chiral 
couplings ci, C3 and 04, the 77 uncertainty dominates the 



errors. Numerical values can be seen at Table Hill Note 
the large discrepancy in the effective range ro for Sets 
I and II with the experimental number. Finite cut-off 
effects are observed in Set III although the 77 induced 
uncertainty would make the value compatible with that 
estimate. Good agreement is observed again for Set IV, 
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TABLE III: Deuteron properties and low energy parameters in the 3 Si — 3 Di channel for the pionless theory (Short), the 
OPE and the TPE potential. We use the non-relativistic relation 7 = ^/2(i np B with B = 2.224575(9). For OPE* we take 
gA = 1-26 as input instead of g^NN = 13.083 as input. The errors quoted in the TPE reflect the uncertainty in the non-potential 
parameters 7, r \ and «o ° m y- Differences from this work are attributed to finite cut-off effects. Experimental values can be 
traced from |54|. 



Set 


Ref. 


7 


V 


As 




Qd 


Pd 


(r- 1 ) 


«o 


a?02 


«2 


ro 


Short 


- 


Input 





0.6806 


1.5265 








00 


4.3177 











OPE 


[41J 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.31(1)% 


0.476(3) 


5.335(1) 


1.673(1) 


6.693(1) 


1.638(1) 


OPE* 


[41] 


Input 


0.02555 


0.8625(2) 


1.9234(5) 


0.2667(1) 


7.14(1)% 


0.471(3) 


5.308(1) 


1.612(1) 


6.325(1) 


1.602(2) 


Set I 


N 2 LO [10J 




0.0245 


0.884 


1.967 


0.262 


6.11 % 


- 


5.420 


- 


- 


1.753 


Set I 


N 2 LO [19] 




0.0256 


0.8846 


1.9756 


0.281 


4.17 % 


- 


5.417 


- 


- 


1.753 


Set I 


This work 


Input 


Input 


0.900(2) 


1.999(4) 


0.284(4) 


6(1)% 


0.424(3) 


Input 


2.3(2) 


3(3) 


1.4(4) 


Set II 




- 


- 












- 








Qo+ TT 


This work 


Input 


Input 


u.yuu^z ) 


i.yyy ^4 ) 




((lj/o 


fi A OA (A \ 


Input 


ooh c\ 
Z.ZZy 10 ) 




1.40(^iy ) 


Set III 


N 2 LO [23] 




0.0256 


0.873 


1.972 


0.272 


5 % 




5.427 






1.731 


Set III 


N 3 LO [29J 




0.0254 


0.882 


1.979 


0.266 


3 % 




5.417 






1.745 


Set III 


This work 


Input 


Input 


0.891(3) 


1.981(5) 


0.279(4) 


7(1)% 


0.436(3) 


Input 


1.88(10) 


5.7(16) 


1.67(8) 


Set IV 


N 3 LO [28] 




0.0256 


0.8843 


1.968 


0.275 


4.51 % 




5.417 






1.752 


Set IV 


This work 


Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1)% 


0.447(5) 


Input 


1.67(4) 


6.6(4) 


1.76(3) 


Nijmll 


[43, 44] 


0.2316 


0.02521 


0.8845(8) 


1.9675 


0.2707 


5.635% 


0.4502 


5.418 


1.647 


6.505 


1.753 


Reid93 


[43., 44] 


0.2316 


0.02514 


0.8845(8) 


1.9686 


0.2703 


5.699% 


0.4515 


5.422 


1.645 


6.453 


1.755 


Exp. 




0.2316 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 






5.419(7) 






1.753(8) 



particularly in the E\ and 3 Di scattering lengths and 
the effective range ro, but only the latter provides a clear 
TPE improvement as compared to OPE. The quantities 
a 02 and 012 are compatible with typical expectations [57j 
from the high quality potential models. 



C. Phase Shifts 

For the a and (3 positive energy scattering states we 
choose the asymptotic normalization 



straints, 



Uk,a(r) 
w k , a {r) 



cosJi — yo(kr) sin^i 



cose A , 
smdi V 

- 1 f ( J2(fcf ) - h{kr) sin <5i) 
smdi V / 



(74) 



1 



sin Si 
tane 
sin<5i 



jo(kr) cos J 2 - y (kr) sin(5 2 



j2{kr) cos 82 — V2{kr) sin^ 



(75) 



where ji (x) = xji (x) and yi (x) = xyi (x) are the reduced 
spherical Bessel functions and Si and 82 are the eigen- 
phases in the 3 S'i and 3 £>i channels, and e is the mix- 
ing angle E±. The use of the superposition principle for 
boundary conditions as well as the orthogonality con- 



u l u 'k,l3 - Uy u k,p + WjW'k,0 ~ w '~i u k,f3 



(76) 



analogous to Eq. (|69|l . to the deuteron wave functions. If 
orthogonality would be applied to the zero energy state 
one obtains an explicit relation of Si, 82 and e with the 
scattering lengths ao, a.2 and 002 as a direct general- 
ization to the coupled channel case the one channel sin- 
glet case given by Eq. H47fl . The explicit expressions are 
rather cumbersome and will not be written down here ex- 
plicitly. The results are depicted in Fig. [5] for Set IV. We 
observe a clear improvement in the threshold region, in 
consonance with the low energy parameters of Table IIIII 
and a moderate improvement over the OPE results in 
the intermediate energy region. This suggests that finite 
cut-off effects may also be built in the phase shifts as well 
as the low energy parameters. 
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FIG. 8: The dependence of the ao2 scattering length (in fm 3 ) 
and the 0:2 scattering length (in fm 5 ) for the TPE potential as 
a function of the asymptotic D/S ratio r\. The boxes represent 
the Reid93 and Nijm II values determined in Ref. |57|l . We 
use Qo = 5.419(7)fm to generate the bands. We use the set 
IV low energy constants. 



VI. ERROR ANALYSIS AND 
DETERMINATION OF CHIRAL COUPLINGS 
FROM LOW ENERGY NN DATA AND THE 
DEUTERON 



A. Propagating experimental errors in ci, C3 and C4. 

The results in the previous sections clearly show that 
deuteron and low energy scattering properties in the 1 Sq 
and 3 Si — 3 D\ channels are sensitive to finite cut-off ef- 
fects and also to the values of the chiral constants after 
removal of the cut-off. We will assume that the values 
for c\, C3 and C4 are free of uncertainties. Then, Set IV 
provides the best description of triplet data but produces 
a slightly off value for the effective range in the singlet 
channel at the 2a confidence level. Note that in the sin- 
glet case the theoretical prediction for r$ does not have 
a large source of error as in the triplet case where un- 
certainties in 77 dominate the error. On the contrary, Set 
III provides a compatible value for the effective range in 
the singlet channel but incompatible values for the triplet 



TABLE IV: Singlet 1 5o and triplet 3 Si - 3 £>i scattering 
and deuteron properties with error estimates using the chiral 
TPE potential. We make a MonteCarlo calculation of the 
input parameters g n NN = 13.1 ± 0.1, Qo, s = —23.77 ± 0.05, 
a 0>t = 5.419 ± 0.007, r\ = 0.0256 ± 0.0004 and the chiral 
constants ci, C3 and C4. The quoted values span an interval 
where 68% of the output is contained. 





Set I 


Set II 


Set IV 


Exp. 


Ci 


-0.81(15) 


-0.76(7) 


-0.81 




C3 


-4.69(1.34) 


-5.08(24) 


-3.20(16) 




C4 


3.40(4) 


4.78(10) 


5.40(1.65) 




ro.s 


9 QQ+0.08 

z - yz -0.04 


9 Q7 +0.03 
*•»<_ 0.02 


2 -°°-0.03 


2.77 ±0.05 


ro,t 


1 S6+ 033 

±.DO_ 75 


1 /I8+ ' 14 
l- 4S -0.25 


1 76+ 003 


1.753 ± 0.008 


A 3 


o.899±S:°8S 


0.900t°;^ 




0.8849 ± 0.0009 


Q d 


0.284l° ;»? 


0-284±g;ggS 


2 76+ 004 

u - z '"-0.004 


0.2859 ± 0.0003 


' m 


1.998™ 


1 QQ8+ 007 
1.9»»_0.007 


l-965±g;gli 


1.971 ± 0.006 


p d 


6.6i8.9 


7 1 +0.9 
'• 1 -0.9 


8.3ll:S 




a?02 




9 9^+0.23 
^• zu -0.16 


1 67+ 013 
— Q.13 






e+3.3 
J - D -6.9 


4 +1<3 

^• u -2.9 


6 71 +0 ' 48 







channel in As and ro at the 3 — 4cr confidence level on 
the experimental side. Thus, on this basis we may reject 
Set III and accept Set IV. 

To improve on this analysis, let us try to include some 
errors on the chiral coefficients. The irN analysis of 
Ref. H3 (Set I) and the NN fit of Ref. [13 (Set II) 
yield some errors. Ref. (Set IV) does not quote 
errors but we will take the educated guess of a 5% er- 
ror for C3 and a 30% error for c\ [5^|. We can propa- 
gate then by a Monte-Carlo simulation implementing also 
the errors in g^NN — 13.1 ± 0.1, ao, s = — 23.77 ± 0.05, 
a ,t = 5.419 ±0.007, 77 = 0.0256 ±0.0004. We assume for 
simplicity that all these quantities are fully uncorrelated. 
This will in general enhance the errors, as compared to 
the case where correlations in c\, C3 and C4 with ttN 
would be taken into account. Perhaps, the best thing 
would be to consider a simultaneous analysis of both 
NN and irN low energy data to build in correlations. 
Obviously, we do not expect good central values for the 
observables judging from Table ITTT1 But there is still the 
possibility of large error bars. 

The outcoming distributions in the low energy and 
deuteron parameters are somewhat asymmetric. Actu- 
ally, for a given set of c\, C3 and C4 distributions we ob- 
serve the appearance of upper bounds in the 3 «Si effec- 
tive range, namely ro,t < 1-79, 1, 75, 1.81fm for Sets I, II 
and IV respectively where the out-coming distributions 
become more dense. The results of the error propaga- 
tion are summarized in Table llVl Thus, we see that the 
values of the chiral coefficients deduced from low energy 
ttN are globally inconsistent, at the la level, with 
the low energy NN threshold parameters after uncertain- 
ties are taken into account. The same remark applies 
to Set II 11]. Again, the loss of predictive power be- 



21 




FIG. 9: Renormalized Eigen Phase shifts for the OPE and TPE potentials as a function of the CM np momentum in the triplet 
3 Si — 3 Di channel compared to the Nijmegen results |43 |. The regular scattering wave functions are orthogonal to the regular 
deuteron bound state wave functions constructed from the OPE with 7 = 0.231605frn -1 , m = 138.03MeV and g-^NN = 13.083 
for the OPE contribution to the TPE potential and g = 1.26 for the TPE contribution to the TPE potential. We take Set IV 
(see main text). 



comes manifest for all the sets although Set IV provides 
the best central values and the smallest errors. The sit- 
uation for the quadrupole moment is noteworthy since 
the difference to the potential value is attributed to Me- 
son Exchange Currents (MEC) and relativistic effects, 
which provide a correction of about O.Olfm 2 (see Ref.[25] 
in Ref. [l^ and also Ref. |59| ). As we see, this is about the 
size of the error deduced from our analysis. In the case 
of the deuteron matter radius the situation is even worse 
since MEC's contributions are much smaller 0.003fm [5^ 
while our predicted errors are larger. It would be ex- 
tremely interesting to reanalyze the problem with the 
present deuteron wave functions |60|| . 



B. Determination of ci, C3 and C4. 

Another possibility is to attempt a direct fit to the 
data. The standard approach is to fit the partial waves 
to a NN database |43L |44| . The problem with such an 
approach is that, unfortunately, there is no error assign- 
ment on the phase shifts and hence a reliable assessment 
of errors cannot be made. Actually, besides the work of 
Ref. 0, where a full partial wave analysis was under- 
taken, other works [2^, H|| |2^ assume fixed values for c\ 
, C3 and C4 without attempting any error analysis based 
on input uncertainties. But even if data for the phase 
shifts with errors were known one expects the quality of 
the fit to worsen as the energy is increased as we think 
that the chiral approach to NN interaction should work 
best at low energies. If the data were known with uni- 
form uncertainty one would fit until \ 2 /DOF exceeds 
one providing an energy window. In such a fit all points 
are equally weighted while we know that the description 
at low energies, where the theory works best, will be com- 
promised by the highest possible energy within such an 



energy window. 

Along the previous line of reasoning we propose, in- 
stead, to fit directly the low energy threshold parameters 
which central values and errors are well known and widely 
accepted. The basic ingredient is to use purely hadronic 
information in the process to avoid any contamination 
due to electromagnetic effects. Specifically, we make a 
Monte-Carlo sampling of the input data assuming that 
as primary data they are gaussian distributed and uncor- 
rected. For any of the samples we make a \ 2 fit to the 
values ro ;S , tqj and As, i.e. we minimize 

Y 2 = ( ro,-r~? \ 2 / r , t -rg? \ 2 / A s - A^ V 
X \ Ar Q , s J + \ Aro. t J + \ AA S J 

(77) 

and determine then the optimal values of c\ , C3 and C4. 
We only accept values where x 2 < lj an d the resulting 
distribution of chiral constants c\, C3 and C4 is given in 
Figs. El As we see, there is a very strong, almost linear, 
correlation between C3 and C4. This can be easily un- 
derstood in terms of the short distance dominance of the 
singlet effective range, since for the pure Van der Waals 
contribution, and in the limit of large scattering length 
r , s ~ R = (MCq) 1 / 4 , with C* 6 given in Eq. JSJ|. Devia- 
tions from linearity are induced from the larger relative 
error of ro, s (1%) as compared to ro,t an d As (0.1%). 
This is different from the large scale partial phase-shift 
analysis obtained in Ref. |2l| where a very small correla- 
tion between C3 and C4 of about 0.2 was found. We have 
checked that cutting-off data with decreasing values of 
X 2 , excludes the points where the distribution is sparse, 
so that the dense part indeed reflects the uncertainties 
in the input data. The fact that the three coefficients 
seem to be on a line is just a consequence of solving by 
minimization a system of three equations and three un- 
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knowns. 

We use g^NN = 13.083, a . s = -23.77 ± 0.05, a 0<t = 
5.419 ± 0.007, r\ = 0.0256 ± 0.0004 and fit a , c 3 and c 4 
to the values r , s = 2.77 ± 0.05, r , t = 1.753 ± 0.08 and 
A s = 0.8846 ±0.0009. Our final result for a sample with 
125 points with x 2 < 1 is 

ci = -1.13t°;^(stat)GeV _1 , 
c 3 = -2.60loi|(stat) GeV -1 , 
c 4 = +3.40lo.4o(stat) GeV" 1 . (78) 

The central value is the mean and the errors have been 
obtained by the standard method of excluding the 16% 
left and right extreme values of the variables, so as to 
have 68% confidence level between the upper and lower 
values. Cutting-off data with % 2 < 0.5 does not change 
significantly the result. 

At the 2(7 level, our values for ci, c 3 and C4 are com- 
patible with the analysis of low energy irN scattering 
of Ref. H3, ci = -0.81 ± 0.15, c 3 = -4.69 ± 1.34 and 
c 4 = 3.40 ± 0.04, but incompatible with the NN full par- 
tial wave analyses j^, where an energy dependent 
boundary condition at a = 1.4fm was used. It is diffi- 
cult to say whether other determinations for the chiral 
couplings based on NN scattering are incompatible with 
ours, since no error estimates have been provided. 



C. Estimate of the systematic errors 

As we have mentioned, any approach based on power 
counting of the potential cannot make an a priori esti- 
mate of the accuracy of the calculation. Nevertheless, we 
can have an idea by simply varying the input parameters. 

At LO we may use either g& = 1-26 as input or 
g-nNN — 13.1, since the difference is the Goldberger- 
Treiman discrepancy, which should be a higher order cor- 
rection. The effect can be seen by comparing OPE with 
OPE* in Table ITTT1 When compared to the TPE result, 
for e.g. Set IV, the error is underestimated this way. 
At NNLO we use the same procedure in the TPE piece. 
Again, the difference should be higher orders. Numeri- 
cally this is equivalent to include g-^NN — 13.1 ± 0.1 and 
consider this systematic error as an statistical one, which 
has already been taken into account. 

We can estimate the systematic error in the chiral con- 
stants by varying the input used to determine the c's. To 
correlate the singlet and triplet channels we must keep 
ro iS and ao. s - So we can interchange the inputs As, ro.t 
with the outputs r m , Qd, cto2 and a-i. This yields a total 
of 15 possible combinations. Another question concerns 
the assessment an error to the fitted variables whenever 
there is no direct experimental quantity, since this choice 
weights the determination of the c's. This is the case for 
«02 and «2, where we make the educated guess of taking 
the difference between the Reid93 and Nijmll values as 
determined in Ref. 1531 as an estimate of the error. The 



situation with Qd is a bit special and we exclude it from 
the analysis 13 . The results are listed in Table [V] We 
see that this estimate on the systematic error provides 
a larger fluctuation than direct propagation of the input 
errors for ci . Symmetrizing the errors we get 

ci = -1.2±0.2(syst)GeV _1 , 
c 3 = -2.6±0.1(syst)GeV _1 , 
c 4 = +3.3±0.1(syst)GeV~ 1 . (79) 

If we attempt a fit to all observables assigning AQd — 
O.Olfm 2 and Aa 02 = 0.4fm we get ci = -0.9, c 3 = -2.71 
and c 4 = 3.85 with a large \ 2 /DOF — 3 basically due to 
the small errors. Obviously a more realistic estimate of 
the errors would be desirable. 



VII. THE ROLE OF CHIRAL VAN DER WAALS 
FORCES 

As we have pointed out, our approach is not the con- 
ventional one of adding short distance counterterms fol- 
lowing a given a priori power counting regardless on 
the approximation where the long distance potential has 
been constructed. Instead, the potential power counting 
dictates the form of the short distance physics by de- 
manding a finite limit when the regulator is removed. In 
order to stress the differences with previous approaches 
it is interesting to see how much of the phase shifts is de- 
termined from the short distance chiral potential without 
adding a short range contribution to the effective range. 
In the standard approach this can be achieved by adding 
a counterterm C2 in the S— wave channels. In Ref. |4l| we 
showed that both perturbatively and non-perturbatively 
the orthogonality constraints for the OPE potential im- 
ply C2 = 0. Here we will see that the bulk of the 5— wave 
interaction can be explained mainly in terms of the chi- 
ral Van der Waals force when renormalization is carried 



TABLE V: Central values for the chiral constants ci, C3 and 
C4 depending on the input, we only include \ 2 < 1- 



Fitted 


Cl 


C3 


c 4 


x 2 


r , s , As,ro,t 


-1.09 


-2.61 


3.36 


0.06 


ro,s,As,r m 


-1.23 


-2.57 


3.34 


0.3 


ro,s,ro,t,r m 


-1.45 


-2.54 


3.26 


0.2 


ro, s , ro,t, ct2 


-1.09 


-2.64 


3.17 


0.4 


ro, 3 ,r m ,a2 


-1.03 


-2.70 


3.26 


0.03 



The discrepancy of potential models to the experimental value 
~ O.Olfm 2 , attributed to MEC's and relativistic effects l59l . is 
about two orders of magnitude larger than the error in the ex- 
perimental number ~ 0.0003fm 2 and the discrepancy between 
potential models ~ 0.0004fm 2 . It is not clear w heth er the dis- 
crepancy can be pinned down with similar errors kill . 
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FIG. 10: Correlation plots for the low energy constants ci, C3 and C4 obtained from the renormalized TPE potential by fitting the 
singlet and triplet effective ranges ro, s = 2.77 ±0.05 and ro,t = 1.753 ±0.08 and the deuteron s-wave asymptotic normalization 
As = 0.8846 ± 0.0009. The dispersion in the data reflects the dispersion in all input parameters, g^NN = 13.11 ± 0.08, 
a , s = -23.77 ± 0.05, a 0> t = 5.419 ± 0.007, 77 = 0.0256 ± 0.0004 which provide \ 2 < 1- 



out, without any additional short distance contribution 
or counterterm. 

For a pure Van der Waals potential of the form 



R 4 

U = —jr, 



(80) 



the zero energy wave function can be analytically com- 
puted |4?J in terms of Bessel functions J v (x) . Normaliz- 
ing to the asymptotic form uo(r) — * 1 — r/aa we get 




The effective range can also be computed analytically j^l 
H from Eq. ^3 yielding 



-4i? 2 



4 J R 3 r(|) 2 , 16RT(I) 2 



3a 
1.39473 R 



3 «q 7T 

1.33333 R 2 

a 



3tt 

0.637318 R 3 

a 2 



(82) 



in agreement with the general low energy theorem of 
Eq. Taking the values of Table[Qfor R = {MC e ) 1/4 

one gets in the singlet 1 S'o channel 



r 0s = 2.39811 - 



3.9418 3.23959 



TO.. 



'0, 



= 2.49192 
= 2.2227- 
= 2.29099 



4.25624 



3.63486 



&0,s 

3.8625 
h 

3.59753 



2.57944 
2.82459 

a Q.s 



(Set I), 
(Set II) , 
(Set III) , 
(Set IV) 



(83) 



The numerical agreement at the ten percent level of the 
«o,s independent term with the full chiral TPE result, 
Eq. (|43() ), is striking 14 . On the other hand, first or- 
der perturbation theory in the OPE potential yields (see 
Sect. A of Ref. [HI) in the form of Eq. JSSJ the result 



'■0.. 



8Mir 1 
1.4369 



5.4789 5.8758 

a 0,s Qq a 



(84) 



Note that the coefficient in 1/<Xq s is slightly better de- 
scribed by the OPE perturbative value than the full OPE 
result (see Eq. Q43jl). a not unreasonable result since this 
coefficient is sensitive to the longest range part of the 
interaction. Likewise, the bulk of the ao-independent co- 
efficient is given just by the most singular contribution to 
the full chiral potential. As we see, for large scattering 
lengths the effective range scales with the Van der Waals 
singlet radius R s = (MCg ijJ 1 ' 4 and not with the pion 
Compton wavelength 1/m, confirming the dominance of 
the short distances singularity in the singlet channel. 

For the triplet channel, the equation cannot be solved 
analytically, and the effective range has a correction due 
to the D-wave (see Eq. t|fT7|l ). Moreover, the scatter- 
ing length is a factor five times smaller than in the sin- 
glet case, so that we do not expect in principle such a 
dramatic agreement. If we neglect the mixing with the 
D-wave and take the R t = (MC^sJ 1 / 4 of Eq. JEJ we 



14 The formula |52J can also be used as a numerical test of the 
integration method and of the numerical solution of the differ- 
ential equations. This is a non-trivial condition due to the rapid 
oscillations of the wave function at the origin. We have checked 
that it is accurately reproduced. 
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FIG. 11: Renormalized Eigen Phase shifts in the 1 So and 3 Si — 3 D\ channel for the pure chiral Van der Waals Ce/r 6 potential 
(Vdw), and the pure NNLO terms compared to the renormalized phase shifts with the same parameters from Table for Set 
IV. We also compare to the Nijmegen database |43j. 
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which, using the triplet scattering length value, ctot — 
5.42 yields r ,t = 1.83, 1.87, 1.75, 1.78 respectively, in re- 
markable agreement with the experimental value. An 
estimate of the mixing effect can be made by using the 
largest van der Waals eigen radius R+ = (MCe.+) 1 ob- 
tained by diagonalizing the interaction at short distances. 
From Table [I] Eq. I|82|l and the experimental value of 
the scattering length we get r = 2.00,2.07, 1.95, 2.05fm 
for Sets I, II, III and IV respectively, accounting for about 
85% percent of the full value. Instead, perturbation the- 
ory for OPE, Eq. @ yields r = 0.62fm, and full OPE 
ro = 1.64. Actually, using the relation 
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we get an explicit correlation between ao, s , ro )S , ao,t and 
ro,t regardless on the numerical values of the chiral con- 
stants C3 and C4 . In the range of physical parameters this 
looks like a linear correlation (see Fig. I12|) between the 
singlet and triplet effective ranges. For r ,t = 1.75 one 
gets r . s = 2.34. 
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FIG. 12: Van der Waals correlation between the singlet and 
triplet effective ranges using the experimental singlet and 
triplet scattering lengths. The point represents the experi- 
mental values. 



To check further the dominance of chiral Van der Waals 
interactions, we plot in Fig.^Jthe phase shifts for a va- 
riety of situations including the pure Van der Waals con- 
tributions, as well as the contribution of the NNLO only, 
which reduces to the previous case at short distances but 
decays exponentially as ~ e -2""- a ^ i on g distances. The 
plots confirm, again, our estimations based on the pure 
Van der Waals potential of the effective range for the 
s— waves, and this is the reason why the triplet s— wave 
is better reproduced than the singlet case for Set IV. Ob- 
viously, by adjusting the effective range changing the chi- 
ral parameters C3 and C4 we could obtain a much better 
description of the data. 

The results of this study show that the singularity of 
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the chiral Van der Waals force is not a feature that should 
be avoided, but instead provides a very simple way to 
describe the scattering data for the s— waves. 

Finally, it is interesting to note, that central waves 
based on taking the chiral limit of the potential are less 
accurately described than the phase-shifts obtained from 
the pure Van der Waals contribution. In this limit the 
singlet 1 Sq channel contains in addition to the Van der 
Waals term a 1/r 5 contribution stemming from the NLO 
TPE contribution. The triplet 3 S 1 ~ 3 D 1 TPE contri- 
bution has a similar structure in addition to the OPE 
tensor 1/r 3 singular short distance contribution. 



VIII. THE TPE POTENTIAL AT NLO: A 
MISSING LINK ? 

In the previous sections we analyzed the renormaliza- 
tion of the NNLO potential. In this Section we analyze 
the NLO in the singlet 1 So and triplet 3 S\ — 3 D\ channels 
and the problem that arises in the latter. We argue that 
similar trends are observed in finite cut-off calculations. 
We also suggest several scenarios on how the problem 
may be overcome. 



A. Convergence in the singlet So channel 



In the singlet 1 St 
tances behaves as 



channel the potential at short dis- 
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The singlet coefficient is negative and, according to the 
discussion in Sect. IIIII one has an undetermined short 
distance phase which can be fixed by using the scattering 
length as input. The effective range in the singlet channel 
is given by 

5.499 



r = 2.122 



(NLO) , (89) 



a a 

which compared with the LO and NNLO results, 
Eq. (|43|) . shows a convergence rate. To show that this 
trend to convergence is not fortuitous we display in Ta- 
ble I VII the threshold parameters of the effective range 
expansion fccot<5 = — 1/ao + rofc 2 /2 + i^fc 4 + + v^k 9, 
depending on the terms kept in the expansion of the po- 
tential given by Eq. Q). As we see there is a clear trend 
to convergence, although the higher order threshold pa- 
rameters display a slower convergence rate since they are 
increasingly sensitive to the shorter range regions. This 
trend is confirmed in Fig. ^] for the phase shift. Ob- 
viously, there is scale separation in the singlet potential, 
and higher order potentials although more singular at the 
origin yield contributions in the right direction. 



B. The problem in the triplet 3 Si — 3 Di channel 



The triplet 3 Si — 3 Di potential at short distances has 
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On the other hand, the diagonalized triplet coefficients 



MC. 



3M(-1 - 10g 2 + 17g 4 
h ~ 256tt 3 / 4 



_ 3M(-l-10g 2 + 47g 4 ) 
MCh - - 256^ ' (92) 

and the mixing angle is given by tan0 = y /2 differ- 
ing by —it as compared to the OPE case 41]. For 
0.5356 < g < 0.8217 one would have an attractive- 
repulsive situation ( see Sect. Illlj) . as in the OPE case 0] 
and in such a case one could take either the deuteron 
binding energy or the 3 S\ scattering length. However, 
for the physical value g — 1.26 one has two short dis- 
tance repulsive eigenchannels, and hence one must take 
the exponentially decaying regular solutions at the ori- 
gin. Let us remind the according to Sect. Mil] finite renor- 
malized results can only be obtained by precisely choosing 



TABLE VI: Convergence of the threshold parameters of 
the effective range expansion fccot<5 = — 1/ao + rok 2 /2 + 
U2 fc 4 + v;ik e ' + Vik a in the singlet 1 So channel depending 
on the successive inclusion of terms in the potential U = 
Ulo + ?7nlo + t^NNLO + • • • • LO means LO alone (and taking 
g-nNN = 13.083 and g A = 1.26), NLO means LO+NLO and 
so on. The only input is the scattering length ao besides the 
potential parameters. For the NNLO case we use Set IV for 
the chiral constants Ci, C3 and C4 given in Table 
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FIG. 13: Renormalized Eigen phase shifts at LO, NLO and NNLO as a function of the CM np momentum k in the singlet 1 So 
and triplet 3 Si — 3 Di channels compared to the Nijmegen results ; 43f for different parameter sets. 



the regular solution at the origin. In this case there are 
no short distance phases, and the scattering lengths, as 
well as the phase shifts are completely determined from 
the potential. The (finite) renormalized results are de- 
picted in Fig. El As we see, the singlet 1 S'o phase-shift 
shows a very reasonable trend, since NLO and improves 
on the LO, and it is improved by the NNLO potential. 
We remind that in the three cases the scattering length 
is exactly the same. However, not completely unexpect- 
edly, the triplet channel results worsen the LO ones. In 
the next subsection we show that if, demanding the stan- 
dard Weinberg counting requires the irregular solution at 
the origin, hence yielding to divergent renormalized re- 
sults. 



C. Finite cut-offs and the Weinberg counting 

The special status of the NLO calculation as compared 
to the LO and NNLO ones has been recognized in previ- 



ous studies in momentum space [lOt where regularization 
was implemented by using a sharp cut-off A. As noted 
by these authors, the allowed cut-off variations at NLO 
are smaller (~ 380 - 600 MeV ) than at LO (~ 700 - 800 
MeV) or NNLO (~ 800 - 1000) but the reasons have not 
been made clear. Let us focus on the triplet 3 Si — 3 D\ 
channel. Within our coordinate space renormalization 
scheme this trend can be easily understood. At LO one 
fixes only one parameter, say ao, and because one has 
attractive and repulsive potentials at short distances the 
system will naturally be driven into the exponential reg- 
ular solution at the origin. Obviously, if one would fix 
some other parameter independently, say ro (or equiva- 
lently using a counterterm C2, and not the one predicted 
by the regular solutions, one would be driven instead to 
the irregular solution, not allowing to remove the cut-off 
in practice. In such a situation one would be forced to 
keep the cut-off finite at the scale where the repulsive 
core sets in. However, at LO the Weinberg power count- 
ing does not allow to fix this additional parameter and 
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one can comfortably reach higher cut-off values. On the 
contrary, at NLO one has two repulsive eigenpotentials 
and one cannot fix any low energy parameter arbitrar- 
ily. Otherwise one would be attracted to the irregular 
solutions at short distances. On the other hand, they are 
attractive at long distances, so that one would expect a 
stability region where the potential becomes flat before 
turning into a repulsive core in both eigenchannels. This 
is exactly what one observes in the NLO calculation of 
Ref. 0] • The occurrence of such a plateau is to some ex- 
tent fortuitous since it is associated to the critical points 
of the potential, and not with some a priori estimate 
on the validity range of the NLO potential. Finally, in 
the NNLO calculation because both eigenpotentials have 
an attractive character one can again increase the cut- 
off since there are no irregular solutions in the problem 
where one can be attracted to. It is very rewarding that 
our coordinate space analysis of short distance singulari- 
ties anticipates when these features can be expected. On 
the other hand this does not imply that finite cut-off cal- 
culations are necessarily wrong, simply that the observed 
features when the cut-off approaches the limit can be un- 
derstood. 

The previous discussion can be illustrated in our ap- 
proach by looking at the short distance cut-off depen- 
dence of the effective range, rg (in fm) and the deuteron 
wave function renormalization, As (in fm -1 / 2 ), in the 
triplet 3 5"i - 3 D 1 channel at LO, NLO and NNLO in the 
standard Weinberg counting as presented in Fig. [21 As 
described in Sect. HIlTl any counterterms can be mapped 
into a given renormalization condition. Once these con- 
ditions are fixed we can ask whether other properties are 
finite or not. Thus, at LO we fix the deuteron binding 
energy (one counterterm), at NLO and NNLO we fix the 
deuteron biding energy, the asymptotic D / S-mixing, 77, 
and the scattering length ao (three counterterms). In all 
cases it is clear that by lowering the cut-off at LO and 
NNLO of the approximation one nicely approaches the 
experimental values. This raises immediately the ques- 
tion whether there is a given value of the cut-off where 
NLO improves over LO. As we see such a region does not 
exist. In addition, although there is a nice and clear trend 
in both LO and NNLO for distances below 0.5fm for As 
down to the origin, this is not so at NLO. So, in this case 
it is not true that low energy properties are independent 
on short distance details, in contrast to the standard EFT 
wisdom. Moreover, Fig. [21 shows explicitly the conflict 
between the Weinberg counting and the remotion of the 
cut-off at NLO because ro and As diverge due to the on- 
set of the irregular solution, as anticipated in our study of 
short distance solutions (see Sect. Illlj) . We have checked 
that this is a general feature on both deuteron and scat- 
tering properties. On the contrary, LO and NNLO have 
a rather smooth limit because in these two cases Wein- 
berg power counting on the short distance counterterms 
turns out to be compatible with the choice of the regular 
solution at the origin. Thus, in the 3 S\— 3 Di channel the 
renormalized solution at NLO in the Weinberg counting 



is divergent while LO and NNLO are convergent. In con- 
clusion, the present analysis shows in a somewhat com- 
plementary manner as done in Sect. IVIII Bl that indeed 
the NLO is problematic, at least non-perturbatively. In 
Sect. IIXI we will see that the problem is not solved if 
the NLO contribution is computed within a perturbative 
framework using the exact OPE-distorted wave basis as 
a lowest order approximation. 



D. The role of relativity and the A resonance in 
the renormalization problem 

The requirement of renormalizability may be regarded 
as a radical step, and renormalized LO calculations de- 
mand violating dimensional power counting on the coun- 
terterms [42j in non central waves such as 3 Po, due to an 
attractive 1 /r 3 singularity. To reach a finite limit the au- 
thors of Ref. 0] propose to promote counterterms which 
in Weinberg's power counting are of higher order. How- 
ever, in their proposal it is intriguing that they choose 
to promote just one counterterm in coupled channels, 
while they could have used a coupled channel countert- 
erm, i.e. three counterterms in total. In the boundary 
condition approach we know from the start how many in- 
dependent parameters must be exactly taken to reach a 
finite and unique limit, the reference to power counting is 
only specified at the level of the potential. Note that the 
power counting in the potential fixes its short attractive- 
repulsive singular character, and this is the origin of the 
conflict of assuming an a priori power counting for the 
counterterms. Finiteness requires that some forbidden 
counterterms must be allowed (promoted) but also 
that some allowed counterterms must be forbidden (de- 
moted). In such a framework, our NLO calculations in 
the 3 <Si — 3 D\ channel lead to finite but nonsensical re- 
sults due to the repulsive-repulsive 1/r 5 singularity (See 
Sect. IVIIIf) . On the other hand, if one fixes as required 
by the power counting the scattering length, the limit 
does not exist because one is driven to the exponentially 
diverging solution at the origin (For instance, Eq. i|(77|) 
gives ro — ► —00). How then can we reconcile finiteness 
with fixing of the parameters ?. As we pointed out al- 
ready, the singular short distance behaviour of the chiral 
potential is in fact a long distance feature which changes 
dramatically when changing the long distance physics. 
Actually, one may reverse the argument and use renor- 
malizability as a selective criterium for admissible long 
distance potentials. In the following we want to provide 
at least two possible scenarios how this might happen, i.e. 
, how modifying the potential at long distances by intro- 
ducing physically relevant information the short distance 
behaviour of the potential changes. 

In the first place, the chiral potential, Eq. (JIJ, was 
derived in the heavy baryon expansion. The short dis- 
tance character may change when not taking such a limit 
since the combination Mr does make the order of lim- 
its ambiguous. A proper treatment of relativistic effects 
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FIG. 14: Cut-off dependence of the effective range, ro (in fm) and the deuteron wave function renormalization, As (in fm -1 / 2 ), 
in the triplet 3 Si — [i Di channel at LO, NLO and NNLO in the Weinberg counting. At LO we fix the deuteron binding energy 
(one counterterm) , at NLO and NNLO we fix the deuteron biding energy, the asymptotic D/S-mixing, r/ and the scattering 
length ao (three counterterms) . We use Set IV of chiral coupling constants. 



requires inclusion of antinucleons in loops, and a satis- 
factory EFT treatment of relativistic effects remains a 
challenging open problem because the standard crossing 
vs. unitarity non-perturbative divorce. On top of this 
one should use a satisfactory relativistic two body equa- 
tion, which necessarily makes the problem fully non-local 
in coordinate space. Nevertheless, there exist "relativis- 
tic" potentials where some of the terms of higher power 
in 1/M than the TPE obtained in heavy-baryon ChPT 
are kept |2(| |53 which have 1/r 7 Van der Waals 
short distance behaviour with attractive-repulsive eigen 
potentials 63] meaning that as in the OPE case one has 
one free parameter. A calculation using these incomplete 
"relativistic" potentials will be presented elsewhere |r3^ |. 

A second scenario is related to the role played by the 
A resonance 15 not included in the present analysis. As 
pointed out in Ref . || , the A provides the bulk of the 
chiral constants, yielding — C3 = 2c4 = g\/2A, with 
A = 293MeV the nucleon-delta mass splitting, yielding 
c 3 = -2.7G6V 1 and c 4 = 1.35GeV~\ The difference to 
the parameters of Ref. |2(j may be due to some other res- 
onances. On the other hand, in terms of scales one has 
A ~ 2m v , which might be regarded as a small parame- 
ter. This obviously does not mean that A vanishes in the 
chiral limit. In the standard chiral counting of the poten- 
tial, Eq. QJ, the combinations c% = Mc±, C3 = Mc^ and 
C4 = MC4 are considered to be zeroth order, but accord- 
ing to the previous argument they could be regarded to 
be enhanced by one negative power. Thus the nominally 
NNLO terms containing C3 and C4 might become NLO 
contributions, and hence changing the repulsive-repulsive 
1/r 5 singularity into an attractive-attractive 1/r 6 one. 
On the other hand, the C3 and C4 contributions of the 



standard NNLO dominate the short distance Van der 
Waals contributions. Actually, much of the NNLO poten- 
tial is built from these terms all over the range. Accord- 
ing to this reasoning our NNLO calculation may be closer 
to a NLO one where the NA splitting is regarded as small 
parameter. In fact, taking NLO+A with — C3 = 2c\4 = 
g 2 A /2A and 77 = 0.0256 one gets A s = 0.8869fm~ 1/2 , 
Q D = 0.2762fm 2 , r m = 1.9726fm and P d = 0.06 in over- 
all agreement with Table ITTT1 It would be rather interest- 
ing to look for further consequences of this A-counting 
at higher orders. The importance of the A in the NN 
problem has been stressed in several works already on 
power counting grounds 0, 01 E3, but the crucial role 
played on the renormalization problem, i.e., the fact that 
the cut-off can be completely removed has not been rec- 
ognized. Our discussion suggests that the momentum 
space cut-off could also be confortably removed in this 
A— counting, unlike the delta-less NLO. 

The two possible scenarios outlined above do not prove 
that the requirement of renormalizability is necessarily 
right, but suggest that looking into the short distance 
singular behaviour of long distance chiral potentials to- 
gether with the mathematical requirement of finiteness 
may provide a significant physical insight into the NN 
problem. In the language of Ref. |42| where promotion 
of counterterms on the basis of the renormalizability re- 
quirement has been stressed, we are perhaps led also to 
the demotion of counterterms (like for relativistic poten- 
tials), or alternatively the promotion of terms in the po- 
tential (like in the A counting described above). 



E. Van der Waals forces, the molecular analogy 
and the chiral quark model 



We thank D. Phillips for drawing our attention to this point. 



The previous arguments show that it is possible to 
change the attractive/repulsive character of the poten- 
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tial at short distances by organizing the calculation of the 
potential in a different manner, but does not give a clue 
on why this actually happens. Remarkably, the analogy 
with atomic neutral systems subjected to Van der Waals 
forces illustrated in Sect. IVfH goes further and provides 
valuable insight into the problem. In low energy molecu- 
lar physics where one works in a Born-Oppenheimer ap- 
proximation, all atomic constituents, electrons and nu- 
clei interact through the Coulomb force arising from one 
photon exchange. At long distances between distant elec- 
trons the potential is a dipole-dipole interaction 



V dip (R) = e 
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where the sum runs over electrons belonging to different 
atoms. In second order perturbation theory the atom- 
atom energy at a separation distance R reads, 



V AA = (AA\V dip \AA) 

+ y, \(AA\V dip \A*A* 

AA^A'A' 



E AA - Ea* 



(94) 



where \AA) and |^4*y4*) is the electron wave function cor- 
responding to a pair of separated clusters in their atomic 
ground state and excited states respectively. The first or- 
der contribution vanishes for atoms with no permanent 
dipole moment. The mutual electric polarization causes 
the Van der Waals interaction between the two atoms, 
Cq /R 6 and because it is second order perturbation theory 
it is obvious that the Cg contribution to the potential will 
always be attractive. However, it is not clear that higher 
order terms would always be attractive. It is remarkable 
that the theorem of Thirring and Lieb [6(j establishes 
that the Coulomb force between constituents implies all 
terms in the expansion being attractive, without appeal- 
ing to the dipole-approximation. Thus, according to this 
result the long distance force will always be singular and 
attractive at short distances, and that is exactly what 
one needs. In such a situation making a long distance 
expansion of the potential, U = — Rg/r 6 — i?g/r 8 + . . . 
and computing the scattering phase shifts by fixing al- 
ways the same scattering length, along the lines pursued 
in this paper, makes much sense. Moreover, one expects 
the results for the phase shifts to be convergent if there 
is scale separation between the corresponding Van der 
Waals radii R$ 3> Rg 3> Our experience with sev- 
eral atomic systems confirms these expectations |67| . 

The argument in the NN system is a straightforward 
generalization of the molecular system above. It is well 
known that there are no colour hidden states between 
colour neutral systems, so that at long distances one may 
assume only exchange of colourless objects. The longest 
range object will be the pion, and the mutual (chiral) 
polarizability will cause attraction between the nucleons, 
exactly in the same way as for atom-atom interactions. If 
we use as an example the chiral quark model, assume for 
simplicity non-relativistic constituent quarks one obtains 



the OPE for quarks. To second order perturbation the- 
ory we get the NN potential in the Born-Oppenheimer 
approximation 



V NN = (NN\V OP e\NN) 

+ y, \(NN\V OP e\HH)\ 2 

HH^NN 
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where Vnn represents the potential in the NN operator 
basis. This yields exactly when HH = NA the results 
found in Ref. j^EllH] and naturally explains why the con- 
tribution from one A intermediate state is attractive at 
short distances. Although this analogy with molecular 
systems is very suggestive, the generalization to all or- 
ders along the lines of the Lieb-Thirring theorem within 
a QCD context remains at present an optimistic specu- 
lation. 



IX. RENORMALIZED PERTURBATION 
THEORY VERSUS NON-INTEGER POWER 
COUNTING 

A. Perturbations on boundary conditions 

In all our calculations we have taken a long distance 
potential calculated perturbatively, and scattering am- 
plitudes have been computed non-perturbatively by fully 
iterating a potential computed in perturbation theory, 
as initially suggested by Weinberg pj. We will call this 
form of solution non-perturbative for brevity. This re- 
quires a non-perturbative treatment of the renormaliza- 
tion problem, which naturally implies that the short dis- 
tance renormalization conditions (or counterterms) are 
determined by the most singular contribution of the long 
distance potential at the origin. For the NN chiral po- 
tential it turns out that the higher the order the more 
singular the potential. As a consequence we have seen in 
Sect. lVIII Cl that for instance Weinberg counting at NLO 
in the 3 Si — 3 D\ channel is incompatible with renormal- 
ization and finiteness due to the short distance repul- 
sive character of the NLO potential. Although naively 
this looks counterintuitive, it is important to realize that 
there are also situations, like LO and NNLO, where 
the regularity condition of the wave function conspires 
against the singularity so that the net effect is well be- 
haved in the scattering amplitudes and deuteron prop- 
erties. Our results in Sect. |^ and \V\ suggest that the 
pattern obtained when comparing LO and NNLO looks 
quite converging numerically, although there appears to 
be no way of making an a priori estimate of the correc- 
tions. 

Perturbative treatments might circumvent this diffi- 
culty since they have the indubitable benefit of allow- 
ing an a priori estimate of the systematic error via di- 
mensional power counting. This causes no problem in 
the calculation of the long distance potential. How- 
ever, we anticipate already that singular potentials are 
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indeed singular perturbations, and power counting may 
not work as one naively expects for the full amplitudes. 
Kaplan, Savage and Wise JH3,|H^| suggested such a per- 
turbative scheme some years ago, where the lowest or- 
der approximation was a contact theory and OPE and 
higher order corrections could be computed in perturba- 
tion theory. This is equivalent to consider mM/ f 2 to 
be first order and m 2 // 2 second order, so that a cal- 
culation involving the chiral constants would be N 3 LO 
in that counting. Unfortunately, the expansion turned 
out to be non-converging at NNLO |70|. In our coor- 
dinate space formulation, this approach corresponds to 
assume for the S waves a boundary condition fixing the 
scattering length ao H3 (See Appendix A of that work) 
and making long distance potential perturbations. In 
our previous work we verified that perturbation theory 
could only account for a contribution to the deuteron 
and 3 5*i — 3 D\ scattering observables at first order. Un- 
fortunately, the second order was divergent, while non- 
perturbatively, i.e., exactly solving the Schrodinger equa- 
tion for the OPE potential, the results were not only fi- 
nite but also numerically quite close to experiment. This 
deserves some explanation. In Fig. 1151 we show the re- 
sults in the deuteron channel when we scale the OPE 
potential C/ope — > A Cope for the s-wave function nor- 
malization As{\) and the effective range fo(A) as a func- 
tion of the scaling parameter A by keeping the deuteron 
binding energy fixed to its experimental value. The non- 
perturbative result is compared to the first order per- 
turbation theory used in Ref. |4l| . Clearly, perturbation 
theory fails even for weak coupling. The experimental 
value for r$ could be obtained by adding a counterterm 
C*2 as done by Kaplan, Savage and Wise [H, |(3jj. A 
non-vanishing C*2 not only violates the orthogonality of 
the zero energy and deuteron wave functions for a long 
distance local potential but also introduces a new pa- 
rameter, reducing the predictive power. Moreover, the 
non-perturbative inclusion of this Ci counterterm with 
the OPE potential yields divergent results (see the dis- 
cussion in Sect. IVIII Cll . In fact, much of the strength 
of C2 is naturally provided by the short distance 1/r 3 
singularity of the OPE potential. 



B. Perturbations on the OPE potential 

Recently, NoggajTimmermans and van Kolck (NTvK) 
have suggested [42| treating the OPE effects non- 
perturbatively, i.e. to all orders, while TPE and higher 
as well as A contributions should be computed in per- 
turbation theory (see also Refs. [55l ITU l72| for related 
ideas). In this section, we analyze such a proposal disre- 
grading the Delta. Our main conclusions will not change 
although numbers could be modified. In such a situation 
the perturbative expansion is equivalent to consider the 
toM/(47t/ 2 ) to be zeroth order while m 2 /(47r/) 2 is taken 
to be second order and m/M is first order, so that the 
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FIG. 15: Dependence of the s-wave function normalization 
As(\) in the deuteron and the effective range ?"o(A) in the 
3 Si — 3 Di channel when one scales the OPE potential (7lo — * 
XUlo- In all cases we fix the deuteron binding energy to its 
experimental value. 
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(96) 



amounts to take g^NN — 13.1 in the OPE piece, hence 
accounting for the Goldberger-Treiman discrepancy. In 
perturbation theory, we must take with gA = 1.26 
and include OPE corrections to higher order. Note that 
the missing first order implies substantial simplifications 
in the perturbative treatment. Indeed NNLO can be 
done within first order perturbation theory since going to 
second order perturbation theory considers which 
is N 3 LO. In the remaining of this section we analyze 
some aspects of such proposal by scaling the strength of 
the perturbation and show that the appearance of non- 
analytical behaviour is intrinsic to singular potentials, 
yielding to perturbative divergences. As we will see, finite 
perturbative calculations including chiral TPE to NNLO 
would require 4 counterterms for the singlet Sq and 6 
counterterms for the triplet 3 Si — 3 Di channel. Our non- 
perturbative results, i.e., fully iterated NNLO potentials, 
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are based on just 1 and 3 counterterms respectively. 



C. Singlet 1 So channel in distorted OPE waves 

Let us examine first the channel, and consider the 
effect of the NLO and NNLO TPE potentials on top of 
the LO OPE potential in long distance perturbation the- 
ory. The fact that they are taken second and third or- 
der respectively means that the effect will be additive at 
NNLO in the scattering properties. For the total poten- 
tial in Eq. (|96[) we write the wave function as 



u k (r) = u^(r)+u^(r) + u^(r) + 



(2), 



,(3), 



and the phase shift becomes 



So 



c(0) 



?(2) 



f(3) 



(97) 



(98) 



At LO the scattering length is a free parameter which 

we fix to the physical value, ag = ao- As we did in our 
non-perturbative treatment in Sect. lVIII Al we will keep 
the scattering length fixed to its experimental value at 
any order of the approximation, so that differences may 
be only attributable to the potential. In the normaliza- 
tion of Eq. (|44() the correction to the phase shift is just 
given by 



s(2) _ 



-A; sin 5, 



U {2) {r)u ( °\r) 2 dr, (99) 



and a similar expression for 5 which can be deduced by 
the standard Lagrange's identity. Here, a short distance 
cut-off r c has been assumed because at short distances 
the NLO potential diverges as C/nlo ~ 1/r 5 . The previ- 
ous formula yields a change also in the scattering length, 
so that we may eliminate the cut-off radius by subtract- 



ing off the zero energy contribution by fixing a. 



(2) 



0. It 



is convenient to recast the result in the form of an effec- 
tive range expansion in the OPE distorted wave basis, 



k cot Sq 



1 

a 



k cot 6, 



(o) 



1 



+ / drU^(r) 4 0) (r) 2 -4 0) (r) 



- / drU^(r) [4 0) (r) 2 -4 0) (r) 2 , 

(100) 



which guarantees = a^' = 0, due to the one sub- 
traction. If we expand in powers of the energy the LO 
wave function we get 



where 



i^"(r) + U(r)u^(r) = 
,(o), 



-u { ° h '(r) + U(r)u { 2>(r) 



2 1 ' 1 ~ ' 1 ' "'_> '• ' ' - u< o\ r ) 

u { ? ) "{r) + U{r)u ( °\r) = u ( ° ] (r) 



(102) 



These equations can be solved recursively. Thus, the 
NLO correction to the effective range is given by 



..(2) 



C/ (2) (r)4 0) W4 0) W^- (103) 



To estimate the short distance contribution we use 
the OPE exchange potential in the form C/lo = 
-e- mr /(R 8 r), with R s = 16f 2 n/g 2 m 2 M the character- 
istic length 1 5o-channel scale. Note that the OPE poten- 
tial in the 1 So channel is Coulomb like at short distances 
for which the complete regular plus irregular solution is 
known. One could then use a short distance expansion of 
the general analytical Coulomb solution. This facilitates 
guessing the solution at short distances for the zeroth 
energy wave function. The higher energy wave functions 
can be computed straightforwardly, yielding for r — ► 
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c r 7 R s + 0{r 9 



(104) 



as can be readily checked by solving the Schrodinger 
equation in powers of energy, Eq. H102(l . The coefficients 
ci and Co correspond to the linearly independent regular 
and irregular solutions respectively and are determined 
by matching to the integrated in asymptotic condition 
u o > \ r ) ~ * 1 ~ r / a o at large distances and at zero energy. 
Obviously, the irregular solution contributes, cq ^= 0, be- 
cause ao is taken to be independent of the potential, and 
hence terms proportional to the coefficient ci are sub- 
leading. Thus, we get 



„(2) 



dr^< ~ 2 



{-4R s )r. 



(105) 



Thus, we conclude that the first order perturbative re- 
sult is badly divergent. This is very puzzling since the 
non-perturbative calculation in Sect. IVHI Al yields a fi- 
nite number (see Eq. and suggests non-analytical 
dependence on the coupling constant. To enlighten the 
situation, let us scale the NLO potential by a factor 
A, E/nlo - * AC/nlOj an d compute non-perturbatively 
the effective range as a function of the scaling param- 



u< k\ r ) = u o°''( r ) + k 2 u^\r) + k 4 uf\r) + . . . . (101) eters, ro(A), with the obvious conditions ro(0) 



,(0) 



and 
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Then one would get 

kcotS = k cot <y<°) + \(ro- r { ^\ Ir 



dr 



uW(r) + uW(r) 



Note that this equation requires assuming ro 



0(A), while non-perturbatively we find r — r ' — 
Now the NLO and NNLO corrections to the vi parameter 
would come as a prediction, 



FIG. 16: Dependence of the effective range J*o(A) in the So 
channel when the scaled potentials U = Ulo + A[/nlo and 
U = Ulo + A(?7nlo + (7nnlo) are considered. In all cases we 
fix the scattering length to its experimental value. 
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2u ( °\r)u^'(r)+u^{r) 



dr 



U i2) {r) + U (3) (r) 
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(107) 



which is also divergent since the leading behaviour of the 
integrand is ~ 1/r 3 at NLO and ~ 1/r 4 at NNLO for 
small r, see Eq. (|104ll . Thus, a further subtraction would 
be needed, predicting the correction to u 3 , 



,( 2 ) 



,(3) _ 



J"o(l) = r^ L °. The result is presented in Fig. ^j] The 
infinite slope at the origin can be clearly seen. Numer- 
ically, we find that for small A -C 0.1, the correction to 
the effective range behaves as tq — Tq ~ \/A, whereas for 
A ~ 1, it behaves as r — 7q ~ A^. This fractional power 
counting A Q , with < a < 1 is evident from the univer- 
sal low energy theorem, Eq. ijBTj)) . for a potential with a 
single scale, U(r) = F(r/R)/R 2 , and appears when the 
short distance regulator is removed 16 . An explicit exam- 
ple is provided by Eq. H82[l when a pure Van der Waals 
potential acts as a perturbation to a boundary condition 
( a contact theory with only ao), since the strength of 
the potential is XMCq — R 4 , but ro contains R ~ A 1 / 4 , 
R 2 — A 1 / 2 and R 3 — A 3 / 4 . It would be interesting to 
predict a priori this non-perturbative non-integer power 
counting analytically for potentials with multiple scales 
as we have done here numerically [6?" 
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Obviously, to prevent the perturbative divergence one 
could subtract an energy dependent contribution, and 
provide the effective range as an input parameter 18 . 



16 In Ref. l4al a potential well was used as a short distance regulator 
which was not removed, and hence the non-analiticity was not 
seen. 

17 Fractional power counting has also been reported to occur also 
in the EFT an alysis of the three body problem for the pionlcss 
theory [siF^l. 

18 This is equivalent to use a short distance energy dependent 
boundary condition on the solution and hence to violate the or- 
thogonality conditions discussed in Sect. ITTT1 



dr 



U^(r) + U i3 \r) 



(108) 



which is logarithmically divergent due to Eq. 1)104(1 . Fi- 
nally, if a fourth subtraction is implemented a convergent 
prediction is obtained for V4 at NLO and NNLO, 



(2) 



„(3) 



dr 



U (2) (r) + U i3 \r) 



2uf\r)ur(r) + 2u^(r)ur(r)+ur(r) 

(109) 

These four subtractions, needed to make a renormalized 
perturbative prediction of the 1 Sq phase shift at NNLO , 
actually correspond to having 4 counterterms, i.e. fixing 
a o, r o,^2 and V3. This result disagrees with the stan- 
dard Weinberg counting (two counterterms at NLO and 
NNLO in the 1 S'o channel). Moreover, besides the loss 
of predictive power as compared to the non-perturbative 
result where only one counterterm is needed the deduced 
renormalized value for V4 is worsened in perturbation 

theory, since vf ' 1 — — 50.74fm 7 , = — 10.45fm 7 and 

(3) 7 
v\ = — 2.88fm . The situation is summarized in Ta- 
ble IVIII where we show our numerical results obtained 
in perturbation theory as explained above and they are 
compared to the Nijmll and Reid93 potential model cal- 
culations (See e.g. Ref. Although these are not 
directly experimental data it is noteworthy that they dif- 
fer by a few percent while the perturbative calculation is 
about a factor 3 larger. The integrals for V4 are rather 



,(°) 



(0), 
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TABLE VII: Threshold parameters of the effective range ex- 
pansion fccot<5 = —1/ao + rok 2 /2 + V2k 4 + v^k 6 + ^fc 8 in 
the singlet 1 So channel in OPE distorted waves perturbation 
theory. We take U (0) = U[°J with g A = 1.26. For the NNLO 
case we use Set IV for the chiral constants ci , C3 and C4 given 
in Table □ 



l S 


LO 


NLO pcrt 


N 2 LO pcrt 


Exp. 


Nijm II 


Reid93 


ao(fm) 


Input 


Input 


Input 


-23.74(2) 


-23.73 


-23.74 


ro(fm) 


1.383 


Input 


Input 


2.77(5) 


2.67 


2.75 


« 2 (fm 3 ) 


-2.053 


Input 


Input 




-0.48 


-0.49 


^3 (fin 5 ) 


9.484 


Input 


Input 




3.96 


3.65 


« 4 (fm 7 ) 


-50.74 


-61.19 


-64.07 




-19.88 


-18.30 



well converging and the matching between the numerical 
solution and the short distance solutions, Eq. (|104f> . is 
quite stable in the region around r ~ O.lfm. 

Thus, in this particular example of the 1 S'o channel 
one sees that our non-perturbative approach based on 
the choice of the regular solutions at the origin predicts 
the phase shift and hence all low energy parameters from 
o>q and the potential as displayed in Table IVfl A per- 
turbative treatment of the amplitude based on OPE dis- 
torted waves requires to fix ao,ro,«2 and W3 at NNLO. 
The phcnomcnological success and converging pattern 
observed when the potential is considered at LO,NLO 
and NNLO is solved non-perturbatively is very encourag- 
ing. The price to pay is to face non-analytical behaviour 
which implies a non-integer power counting. The trend 
observed here can be generalized to other channels. A 
more thorough discussion of this issue will be presented 
elsewhere 671. 



Here, M is the effective range matrix defined by its rela- 
tion to the unitary S-matrix, 

DSD" 1 = (M + i/cD 2 ) (M - i/cD 2 ) 1 , (113) 

and D = diag(l,/c 2 ). The reduced Bessel functions ma- 
trices are given by j = diag(j , h) and y = diag(yo, £2) 
with ji(x) = xji(x) and yi(x) = xyi(x). At low energies, 
one has the effective range expansion(see e.g. [57) and 
references therein), 



M = -(a)" 1 + -rk 2 + vfc 4 



(114) 



Here, we have introduced the scattering length matrix, 



a «02 
a 2 a 2 , 



(115) 



the effective range matrix 



r = 




(116) 



and so on. These parameters have been determined in 
pS7j from the potentials of Ref. Q • Proceeding similarly 
as in the one channel case, one gets, after one subtraction 
at zero energy the effective range function in perturbation 
theory, 



M + (a)" 1 =M(°) + (a^)- 1 
dr 



u (0)t U (2) u (0)_ u (0)t U (2) u (0) 



(117) 



D. Triplet 3 Si 



D\ channel in distorted OPE 

waves 



We turn now to the triplet 3 Si — 3 D\ channel. The 
reasoning is a straightforward, although tedious, coupled 
channel generalization of the So case, with the addi- 
tional feature that the short distance behaviour is dom- 
inated by a 1 /r 3 singularity (instead of 1 jr) , and so the 
short distance behaviour is different. It is convenient to 
introduce the potential matrix as 



U(r) 



_ Us Sl (r) U El (r) 



U El (r) Us Dl (r)J ' 
and the matrix wave function, 

u k , a {r) u ki p(r)\ 



Ufc(r) = A 



W k ,a( r ) w k,/3(r) 



(110) 



(111) 



with A a constant energy dependent matrix, subject to 
a slightly different normalization than Eq. I|75|) . 



Mr) -> T j(fcr)D _1 M - y(kr)T> 



(112) 



The condition a 



(0) 



ao must be imposed, since a§} an d 



c4°' are predicted from (at LO one only needs one 
counterterm). This formula implies that one introduces 
two new conditions to fix now CV02 and a 2 to their ex- 
perimental value. Along similar lines as done before, we 
analyze the finiteness of the previous expression by com- 
puting the effective range matrix. To this end we expand 
the coupled channel wave function in powers of momen- 
tum 



ui 0) (r) = (r) + k 2 uf \r) + k 4 u { ° ] (r) -f 
to get 

roc r 

u( 0)t U( 2 )u( 0) +u( 0)t U( 2 )u : 



J2) 



dr 



(118) 



(119) 



The LO OPE short distance behaviour of the triplet wave 
functions has been worked out in our previous work |4l| . 
It is convenient to define the triplet length scale 



Rt 



3gjM 
32^/ 2 



(120) 



34 



which value R t — 1.07764fm. One has the general struc- 
ture 



u ( r ) = ~Jk [-5- 



1 / r 



3/4 



+4V2J ^ 



w(r) = 



- C lR f lR (r)e 

C 2fl /2fl(r) e -W ^ + V2C 1A f 1A (r)e- 
V2C 2A f 2A (r)e 4 ^ ^" 

V2C 1R9lR (r)e +i ^ 

V2C 2R92R {r)e-^ % 
C 2A g 2A {r)e il ^l ~ 



1 / r 
73 



3/4 



CiA5iA(r)e^ 



4i> 



(121) 



where the constants Ci#, C 2 fl, Cia and C 2 a depend on 
the energy and the OPE potential parameters. The reg- 
ular solution is selected when one takes C\ R = 0. The 
functions appearing in this formula are of the form 

%/2 



m - 5> i 

n— v 

<m-5>(£) 

n— v / 



n/2 



(122) 



For the present calculation we only need the power be- 
haviour (see Appendix B of Ref. |41j ) 



uf(r) 



r 3/4 

r 3/4+5/2 
,,3/4+5 



(123) 



which shows that, again, the first order correction 
to the effective range matrix is logarithmically diver- 
gent because the NLO potential diverges as 1/r 5 and 
Tj( 2 )uoU2 ~ 1/r 19 . As previously, the situation could 
be amended by adding 3 new counterterms to fix the 
effective range matrix r and then v would come as a pre- 
diction. So, at NLO in perturbation theory one needs 
a total of 6 counterterms to generate a coupled chan- 
nel finite amplitude. When adding the NNLO contribu- 
tion this number of counterterms remains the same since 
U^u^uf ~ r 1 ' 2 and U^uf ] 2 ~ r 1 ' 2 . 

An illustration of non-analytical non-perturbative be- 
haviour in the 3 Si — 3 Z?i-channel can be looked up in 



19 There is a subtlety here. The terms containing the regular ex- 
ponential at the origin are convergent, regardless on the power 
of r in the denominator. Naively logarithmically divergent inte- 
grals would become convergent when combined with oscillating 
functions. However, these functions appear squared so that the 
logarithmic divergence prevails. 
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FIG. 17: Dependence of the s-wave function normalization 
As(X) in the deuteron and the effective range ro(A) in the 
3 Si — 3 Di channel when one scales the TPE potential U = 
Ulo + A(C/nlo + I/nnlo). In all cases we fix the deuteron 
binding energy, the asymptotic D/S-ration, r\ and the s-wave 
scattering length, ao are fixed to their experimental values. 



Fig. 1171 There, the behaviour of the s-wave function nor- 
malization As(X) in the deuteron and the effective range 
ro (A) in the 3 Si— 3 Di channel when the TPE potential is 
scaled as U = Ulo + A(Z7nlo + C^nnlo) and the deuteron 
binding energy, the asymptotic D/S-ration, 77 and the s- 
wave scattering length, ao are fixed to their experimental 
values. 



E. The deuteron in distorted OPE waves 

To conclude our analysis of perturbation theory we 
study now the deuteron bound state. According to 
Fig. 1171 there appears some tiny non-analyticity for very 
small couplings in the asymptotic s-wave normalization 
As- Note that there is an apparent linear behaviour with 
the exception of the very small A region, making one sus- 
pect that the result might be obtained in perturbation 
theory. We will see below by an explicit perturbative 
calculation that this is not so. We have checked that 
this trend also occurs for other quantities such as the 
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quadrupole moment, Qd, the matter radius, r m , and the 
D-state probability, Pq. Here, we show, as it has been 
done above for the scattering problem, that this can be 
traced to a first order divergent renormalized result. 
We define the two component deuteron state, 



u 7 (r) 



(124) 



u 7 (r) 

In perturbation theory, we expand the potential 

U(r) = U<°) (r) + U {2) (r) + U< 3 > (r) + . . . (125) 

and thus the deuteron wave function for fixed energy (or 

7), 

u 7 (r) - (r) + uf (r) + u< 3 ) (r) + . . . (126) 
where (it 7 °^(r), w 7 °^(r)) correspond to the lowest or- 

(2) (2) 

der solutions of the problem and (w 7 (r),«; 7 (r)) and 
(w 7 3 ^(r), w^^)) satisfy 

- uf(r) + [U(°>(r) + 7 2 1 uf(r) = 

U(°)(r)+ 7 2 

1 (3) (f)+ u(0) (r)+7 2 



-u$°(r) 



,( 2 ) 



,( 3 ) 



(r) = -U( 2 )(r)u' 0) (r) 



(r) = -U< 3 )(r)uf (r) 
(127) 

We look for normalized solutions, so that perturbatively 
1 = 



druf t( r ) u (o) (r) 

dr(u( 2 )t (r)u (0) (7 . ) + u (0)t (r)u (2) (r) 
" dr (uft(r)uf (r) + uf( r ) u ( 3 )(r) 



(128) 

The zeroth order equation was solved in our previous 
work |41|. where it was shown that 7 was a free param- 
eter, which means 7*- -' = 7, and the regular solution at 
the origin was selected (see Eq. (I121f) l to ensure normal- 
izability at the origin. We will keep always the same 
fixed value at any order of the approximation, so that 
0. To analyze the NLO and NNLO prob- 



7 



(2) 



7 



(3) 



lem analytically we proceed by the variable coefficients 
method. The zeroth order equation is a homogenous lin- 
ear system with four linearly independent solutions, 



,(0) 



(r) 



Ui(r) 
Wi(r) 



1,2,3,4 



(129) 



The first order equation is an inhomogeneous linear sys- 
tem, which solution can be written as 



uf(r) + uf (r) = £>(r)ui(r) 
c l (r)w t {r) 



w^{r)+wf(r) 



,=1 
4 

E< 

i=l 



(130) 



The variable coefficients satisfy 
4 

5>KrK(r) = 
i=i 
4 

^c-(r) Wi (r) = 

i=l 
4 

5><(rK(r) = F u (r) 
1=1 

£cJ(rX(r) - F w (r) (131) 
1=1 

where we have defined the driving term 

P(r) = ( K{r) ) 

= -U( 2 '( r )uf( r )-U' 3 '(r)uf(r) (132) 
which at short distances behaves as 

F(r) - r 3/4 C 5 r 
whereas at large distances one has 



- r^ 4 C 6 r- 6 



(133) 



F(r) ~ e^ r e - mr 



(134) 



To proceed further, we choose the following linearly in- 
dependent solutions fulfilling the asymptotic boundary 
condition at infinity 



ux(r) - 




wi(r) - 


- 0, 


u 2 (r) - 


- 0, 


w 2 (r) - 


V i r 


U3(r) - 




w 3 (r) - 


- 0, 


u 4 (r) - 


- 0, 


iu 4 (r) - 


V 7^ 



( 7r ); 



(135) 

Any of these solutions has a short distance behaviour of 
the general form given in Eq. (|121fl . So that, all these 
solutions are necessarily singular at the origin. Using 
Krammer's rule the solutions to the linear differential 
system, Eq. <|131[) . which are regular at infinity read 



ci(r) 
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W 



dr 



U2 U 3 U4 

W2 W3 W4 

F u u' 2 u' 3 u 4 

F,„ w'o w'» w' A 



(136) 



3G 



i r 

C2{r) = wJ 



dr 



Ui u 3 U4 

wi W3 W4 

u[ F u u' 3 u\ 

w[ K w 3 w 'a 



(137) 



However, the leading short distance behaviour of the in- 
tegrand is given 



cj(r) ~ r 3/4 (e 4V5flA : )2 



r (145) 
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where W is the Wronskian 
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At asymptotically large distances we have 
w (2) (r) 



(141) 



and similarly for the N 2 LO correction. Note that the 
normalization condition, Eq. 1)128(1 . implies a linear rela- 
tion between Cg and as well as cb and c^ 1 . The 
total D/S ratio obtained by including the zeroth order 
contribution is given by 



rj = 7} 



(0) 



1 



(142) 



If we fix rj we get a relation between C5 and C£>. The 
coefficients and are given by 



V 



(0) 



c (2) 



,(3) _ 



,(3) 



ci(oo) 
c 2 (oo) 
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The long distance behaviour of the integrands is well be- 
haved since, up to inessential powers in r, one has 



4(r) 
c' 2 (r) 
c' 3 (r) 
c' 4 (r) 



-2mr 



-2mr 



-(2m+27)r 
-(2m+2y)r 



(144) 



So, we expect the coefficients cs and cd to diverge if the 
short distance cut-off is removed, r c — > 0. It is unclear 
how this divergence might be avoided. Unlike the scatter- 
ing problem in perturbation theory, where energy depen- 
dent (and hence orthogonality violating) subtractions are 
needed, it would be difficult to accept a bound state not 
normalized to unity unless one includes, besides pn other 
Fock state components, such as pnn. This short distance 
analysis holds also when the 1/r 6 A-contributions are 
taken into account perturbatively. 

Thus, perturbation theory on the distorted OPE ba- 
sis for the deuteron makes sense only as a finite cut-off 
theory. In the appendix [S] we develop further such an 
approach to NLO, where r\ is an input and to NNLO, 
where both As and rj should be fixed. We also show that 
in the cut-off theory the NLO yields tiny corrections to 
deuteron properties whereas the NNLO dominates. This 
proves that, at least perturbatively and in the absence of 
the A, the (integer) power counting to NNLO in deuteron 
properties is obviously not convergent. To some extent 
this result resembles qualitatively the findings of Ref. [7(j 
based on the idea that OPE and TPE can be included 
perturbatively 68]. Of course, an amelioration of the 
convergence in the cut-off theory when A's are included is 
not precluded, and deserves further investigation. How- 
ever, the very need of a finite cut-off will still hold, as 
our analytical study shows. 

In Fig. [E|we show LO, NLO and NNLO order wave 
functions when a finite short distance cut-off r c = 0.5fm 
is considered. The strong divergence of the wave function 
at the origin can be clearly seen. 



X. CONCLUSIONS 

In the present work we have extended the coordinate 
space renormalization of central waves in NN interaction 
discussed in our previous work ^lj for the OPE poten- 
tial to the TPE potential. As we have stressed along the 
paper, the main advantage of such a framework is that 
the (renormalized) potential is finite everywhere except 
at the origin where a Van der Waals attractive singular- 
ity takes place. This suggests using a radial cut-off which 
provides a compact support for the short range part of 
the potential thus making scheme dependent contact in- 
teractions innocuous for the long range solution. Thus, 
there is no need to device different regularization meth- 
ods for the potential and the wave functions both for 
bound state and scattering state solutions. As a result 
model independent long range correlations between NN 
observables can be deduced if the potential is iterated to 
all orders. 
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FIG. 18: Deuteron wave functions at LO, NLO and NNLO in 
perturbation theory when a short distance cut-off r c = 0.5fm 
is considered for the perturbative corrections. At LO we fix 
7, at NLO we fix 7 and r\ and at NNLO 7, r\ and As are 
fixed. Solutions are perturbatively normalized. The LO wave 
functions are normalized by taking As = 1 



Important constraints can be deduced from the re- 
quirement of a small wave function in the unknown short 
distance region. As a consequence the boundary con- 
dition for the wave function at short distances becomes 
energy independent if the long range contribution to the 
potential is also energy independent. We stress here that 
such requirements, although quite natural from a physi- 
cal viewpoint, may not appear obvious within the so far 
established EFT framework, and it would be very inter- 
esting to provide further arguments within EFT itself 
supporting our unconventional framework . Actually, 
we find that the singularity structure of the potential at 
short distances determines uniquely how many parame- 
ters must be regarded as unknown, non predictable, in- 
formation. This is done in terms of short distance phases 
or cquivalently via suitable mixed boundary conditions 
at the origin. Moreover, for an energy independent po- 
tential the orthogonality of wave functions precludes a 
possible energy dependence of the boundary conditions. 
In the particular cases studied in this paper, namely 1 So 
and 3 Si — 3 D\ channels, we have found that besides the 
NNLO TPE potential parameters, one can use the S- 
wave scattering lengths in both channels as well as the 
deuteron binding energy and the asymptotic D/S ratio of 
the deuteron wave functions as independent input infor- 
mation. The remaining scattering or bound state proper- 
ties in the triplet channel are then predicted unambigu- 
ously. Based on the superposition principle of boundary 
conditions, we have found analytical and simple univer- 
sal rational relations which clearly exhibit these features. 
These universal relations would be very difficult to de- 
duce in momentum space and, moreover, are free from 
uncertainties attributable to finite cut-off effects. So, the 
cut-off has been effectively eliminated. On a numerical 
level, the fact that our problem is an initial value problem 



for the Schrodinger equation starting at infinity, makes 
possible to obtain any solution by competitive algorithms 
with adaptable integration steps with any prescribed ac- 
curacy. This allows to faithfully describe the short dis- 
tance oscillations of the wave function. This is in contrast 
to the standard Lippmann-Schwinger treatments, where 
matrix inversion methods may eventually run into com- 
puter space limitations with a natural loss of space res- 
olution as a side-effect. The non-trivial oscillating struc- 
ture of the wave functions with ever decreasing periods 
of the wave functions close to the origin would actually 
be very difficult to reproduce within a momentum space 
framework. 

According to our analysis, there are finite cut-off effects 
in previous works dealing also with TPE potentials both 
in coordinate as well as in momentum space. The induced 
corrections are larger than the experimental uncertainty 
of the computed observables, so that in some cases agree- 
ment with data may be clearly attributed to the choice 
of a finite cut-off. In our energy independent boundary 
condition treatment we found short distance cut-offs of 
about r c — 0.1 — 0.2fm to be rather innocuous. Within a 
Wilsonian viewpoint of renormalization, changes in the 
cut-off should correspond to decimation, i.e. halving, and 
not to linear changes in the scale. If one associates this 
coordinate space cut-off to a momentum space ultraviolet 
cut-off of A = 7r/2r c we are dealing with an equiva- 
lent momentum scale of about 1.5 — 3GeV, much larger 
than the scales below lGeV usually employed in momen- 
tum space calculations where only linear sensitivity to 
changes of the cut-off was implemented. Nevertheless, it 
is fair to say that the calculations based on Sets III and 
IV provide not too large discrepancies. 

As one naturally expects in a renormalized theory, er- 
rors are dominated by uncertainties in the input data, 
and not by cut-off uncertainties. Indeed, we seem to 
reach a limit in the accuracy of the predictions, parallel- 
ing the findings in ChPT for mesons at the two loop level. 
At the OPE level, one can predict bound state and scat- 
tering properties in the singlet 1 Sq and triplet 3 Si — 3 D\ 
channels solely from the deuteron energy and the 1 Sa 
scattering length. At the TPE level, one needs not only 
the additional chiral constants c\ , C3 and C4 but also the 
triplet S-wave scattering length and the asymptotic D/S 
deuteron ratio. Although the TPE central values predic- 
tions improve, the induced TPE errors turn out to be 
larger than the OPE uncertainties. In fact, this large un- 
certainties make that, within errors, the TPE calculation 
becomes compatible with experimental data at the la 
level. This suggests that in order to see in a statistically 
significant sense other effects, such as electromagnetic, 
relativistic and three pion effects one must first improve 
on the input data. Otherwise, predictive power is lost. 
Nevertheless, given the finite cut-off effects detected in 
previous works, the role of these corrections beyond TPE 
should be reanalyzed within the present approach. 

One of the important consequences of our treatment is 
that the chiral constants c\ , C3 and C4 can be determined 
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from low energy data and deuteron properties. Specifi- 
cally, we have used the singlet and triplet effective ranges 
as well as the asymptotic S-wave deuteron wave function 
to c\ , C3 and C4 with errors varying all input data within 
their experimental uncertainties. The decision on what 
set of data should be used to pin down the chiral coef- 
ficients is not entirely trivial, because it should become 
clear which hypothesis we want to verify or to refute. 
The absence of cut-off effects makes this test cleaner; we 
just check whether the TPE potential holds from zero 
to infinity. Obviously, this cannot be literally true, but 
one expects that at low energies other short range effects 
can be considered negligible. Let us remind that error 
analysis within NN calculations was only carried out in 
a large scale partial wave analysis to data in Ref. [Tlj . 
The determinations of chiral constants based on a fit to 
NN databases 13 E3 f° r phase shifts lack any error 
estimates because the databases themselves are treated 
as errorless. The determination of chiral constants from 
peripheral waves has similar drawbacks. From the chiral 
theory point of view we see that it is possible to deter- 
mine these parameters precisely in the regime where we 
trust the theory most, namely in the description of low 
energy NN data. A fit becomes possible, and the values it 
yields only differ by 2a with the determination from irN 
data. We do not exclude that our values for the chiral 
constants may eventually spoil the successful overall fit of 
phase shifts in all channels presented in the past, after all 
renormalization has been carried out. If so, the situation 
on the effectiveness of effective field theory would be in a 
less optimistic shape than assumed hitherto. A prelimi- 
nary analysis of the problem shows what Van der Waals 
coefficients in the TPE potential correspond to attractive 
short range interactions, and hence what phase shifts are 
completely determined in terms of coupled channel scat- 
tering lengths. This issue is very relevant and would re- 
quire a detailed channel by channel analysis and renor- 
malization, taking as input the scattering lengths docu- 
mented in our previous work |57| and integrating in from 
large distances along the lines of the present approach. 
Full details are reported elsewhere [74| . 

Nevertheless, despite of the good convergence in the 
1 S'o channel for LO, NLO and NNLO calculations, we 
have noted a difficulty for the triplet 3 Si — 3 Z?i chan- 
nel at NLO of the potential. In contrast to dimensional 
power counting expectations one cannot use the scat- 
tering length «0j the effective range ro and ckq2 as ar- 
bitrary input parameters at NLO in the potential (one 
could equally take 7, 77 and cxq) but they are entirely 
predicted from the potential as required by finiteness of 
the phase-shifts. Otherwise, the scattering amplitude di- 
verges, as we have shown. We have also seen that even if 
one assumes a finite value of the cut-off the NLO is worse 
than the LO, suggesting that the problem may indeed be 
related to the power counting on the long distance po- 
tential. Remarkably, these parameters must be fixed at 
NNLO where, according to the standard approach, no 
further low energy parameters should be fixed. This mis- 



match in orders can be understood if one considers the 
NA splitting to be a small parameter, making much of 
the NNLO contributions to the potential become NLO 
ones, because C3 and C4 would be order minus one. In 
such a case, our interpretation goes hand in hand with 
the standard approach; one needs three independent low 
energy parameters at NLO in this counting. The con- 
sequences of this A-counting to higher orders within the 
context of renormalization will be explored elsewhere. Of 
course, we should point out that despite the rather tan- 
talizing description achieved at NNLO the existence of 
a consistent power counting guaranteeing the success of 
the present approach to all orders remains to be proved. 
A key ingredient of such a power counting would be the 
correct incorporation of all long range physics. Appar- 
ently, within Weinberg's power counting the NLO in the 
deuteron channel misses important contributions. 

Finally, we have analyzed the consequences of a per- 
turbative expansion of TPE effects taking the OPE re- 
sults as a zeroth order approximation as suggested re- 
cently 01 ■ Our non-perturbative calculations based on 
iterating a perturbative potential to all orders exhibit 
unequivocal non-analytic dependence on the expansion 
parameter, due to the singular character of the chiral po- 
tentials at the origin. This is equivalent to a non-integer 
enhancement of the power counting A" with < a < 1 
in the potential, and it would be interesting to know the 
general rules of such a counting a priori |67| . Thus, per- 
turbation theory based on standard power counting be- 
comes divergent and can only yield finite results at the 
expense of introducing more perturbative counterterms 
than are needed in a non-perturbative treatment. This 
is just a manifestation of the fact that singular poten- 
tials require infinite counterterms in perturbation the- 
ory, while only a few ones are needed non-perturbatively. 
Specifically, our analysis shows that it would be neces- 
sary to include at least 4 counterterms for the singlet 
1 5o and 6 counterterms for the triplet 3 Si — 3 D\ chan- 
nel at NNLO. This proliferation of counterterms is ex- 
pected to occur also in other partial waves because the 
singularity of the potential dominates over the centrifugal 
barrier at short distances. In the 1 Sq channel, we have 
seen that adding more counterterms in fact worsens the 
results for the effective range expansion parameters. In 
contrast, our non-perturbative calculations are based on 
just 1 and 3 counterterms respectively. The good quality 
of our results suggests that our choice of less countert- 
erms cannot be refuted on the basis of phenomenology. In 
the deuteron case we have made a calculation to NNLO 
in perturbation theory. Our analysis shows that such a 
perturbative approach only makes sense if a finite cut-off 
is introduced. In any case, the cut-off theory has less pre- 
dictive power, does not provide a better phenomenologi- 
cal description of the deuteron than our non-perturbative 
renormalized results and is non-convergent since NNLO 
corrections are numerically much larger (two or three or- 
ders of magnitude) than NLO ones, despite being para- 
metrically small. In our view this is a perturbative mani- 
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festation of the short distance dominance which has been 
unveiled non-perturbatively. In addition, the difficulties 
faced by a perturbative treatment are simply absent in 
the non-perturbative approach. 

One of the main goals of nuclear physics is the de- 
termination of the nucleon-nucleon interaction. From 
a theoretical viewpoint the disentanglement of such an 
interaction in terms of pion exchanges based on chiral 
symmetry requires dealing with non-trivial and, to some 
extent, unconventional non-perturbative renormalization 
issues in the continuum, but it is crucial because it shows 
our quantitative understanding of the underlying theory 
of quarks and gluons in the chirally symmetric broken 
phase. Our results also show that the singular chiral Van 
der Waals forces are not necessarily spurious and inconve- 
nient features of the chiral potential. Instead, as we have 
shown, the singularities alone in conjunction with renor- 
malization ideas explain much of the observed S— waves 
phase shifts with natural values of the chiral constants, 
and provide an appealing physical picture. In this regard, 
it is interesting to realize that based on the analogy with 
molecular systems, which also exhibit a long range Van 
der Waals force, the liquid drop model was formulated 
more than 60 years ago. Chiral dynamics may provide 
not only a closer analogy and perhaps more quantitative 
insights into the hydrodynamical and thermodynamical 
properties of nuclei but also a theoretical justification 
from the underlying theory of strong interactions. 
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APPENDIX A: THE DEUTERON IN 
OPE-DISTORTED PERTURBATION THEORY 
WITH A CUT-OFF TO NNLO 

In this appendix we illustrate the situation discussed 
in Sect. lIXl bv solving numerically the set of perturbative 
equations H127fl . As we have mentioned such a calculation 
makes only sense within a finite cut-off scheme. 

In practice, we integrate from large distances (~ 25fm) 
with the conditions specified by Eq. (|141|l with some pre- 
scribed values of eg and cd 20 ■ This can be advanta- 



geously done using the superposition principle of bound- 
ary conditions, Eq. (|6U[) . yielding in perturbation theory 

u 7 (r) = u^(r)+u^(r)+u^{r) + ... 

«, 7 (r) = u^(r)+u^(r)+u^(r) + ... (Al) 

At LO the wave function can be written as 

«4°>(r) = 4 0) (r)+7 ? W U g ) (r) 

w W(r) = wPi^+ri^w^ir) (A2) 



and is determined from the regularity condition at 
the origin plj . At LO the normalization factor is 



dr(u^(r) 2 +w["\r) 2 ) (A3) 



(0)\2 



The NLO and NNLO contributions are 

u P0(r) = C fuf{r)+^c^u^\r) 

uf{r) = cfuf(r)+r,^c^u%\r) 

w ( 3 \r) = cfw^i^+rj^c^w^ir) (A4) 



The advantage is that the functions appearing here only 
depend on the potential and the deuteron binding energy, 
whereas the coefficients must be determined by some ad- 
ditional conditions. In the first place, normalization to 
NLO and NNLO requires orthogonality of the wave func- 
tions to the LO solution, 

/>oo 

= / dr(J°)(rV 2 >(r) + k/ 0) (r> (2) 0)) 



0=1 dr(u 

(0) ( r ) M (3) ( r ) + w {0) ( r ) w (3) ^ ( A5 ) 



This implies the couple of linear relations 21 



,(0) 



£ { d _ JiT rfr(u(°) (r)4 2) ( r ) + ^ (0) ( r ) w s } ( r )) 
4 2) ~ / O ° dr (« (0) (» , )«D ) ( r )+«' (0) (» , )wD ) (r)) 

c d ] io°° dr ( ui0) ( r )4 3) ( r ) + w(a) ( r ) w s ] M) 

(A6) 



Further relations can be obtained by imposing renormal- 
ization conditions. Note that the required number of con- 
ditions increases with the order. This is similar in spirit 
to the procedure of adding more counterterms for the 



direct use of the explicit expressions Eq. 11431 involving deter- 
minants. 

21 For instance, at r c = 0.5fm we get cS' = — 5.865cg 2 ' and Cj^ = 



This is numerically more efficient and stable procedure than a 



-5.545c 



(3) 
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scattering problem discussed in Sect. IIXI For instance, 
using the perturbative expansion for Ag and Ap 



for r c > 0.2fm. At NNLO we fix 77 and As- The solution 



As = 4 0) ( 



1 



..(2) 



„( 3 ) 



-4, 



= 4°V 0) (l + cg> + cg> +...)=„* 



(A7) 



In practice, we use a short distance cut-off r c for the NLO 
and NNLO contributions only. Deuteron properties can 
be written to NNLO as follows 

r - r (0) + c (2) (2,S) , (0) (2) (2,£>) 

+ c?V^+^cgV^ + ... (A8) 

- ^ 0) +4 2) q^ s) +^4 2 )q^ 

+ cfQ^+^c^Qf D ^ + ... (A9) 

where the potential contributions have explicitly been 
factored out. The numerical solution requires some care, 
due to the short distance instabilities and oscillations. 
This requires using an adaptive grid to optimize the con- 
vergence. Since solutions of different orders must be 
mixed in the evaluation of the orthogonality conditions, 
Eq. HA6fl . and observables, Eq. I|A9[) . we solve all LO, 
NLO and NNLO equations simultaneously to provide all 
functions on the same grid. 

At NLO and fixing r = r c we demand the experimen- 
tal value of r], from Eq. (|142|l and Eq. (|A6(I . This way, a 

(2) (2) 

solution which we denote by (c s |nlo, c~d |nlo) can be 
obtained. From there we can obtain deuteron properties 
to NLO, as a function of the perturbative cut-off r c In 
Fig. E| we show the dependence of As, r m , Qd and pd 
on r c . As we see the NLO correction is tiny and stable 



is now c 



(2) 



(2) 



S iNNLOjCp I NNLO 



and 



,(3) 



(3) 



S InNLOjCq I NNLO 



(2) (2) 

Note that in general the NLO coefficients c s and c D 
must be readjusted. In this case the correction is much 
larger than the NLO case (see Fig. I19J) and the cut-off 
dependence is stronger due to the the 1/r 6 singularity of 
the NNLO potential. As a curiosity, we mention that at 
short distances worry-some negative D-wave probabilities 
show up below r c = 0.1 7fm at NNLO, a spurious feature 
which can only take place in perturbation theory and 
sets a unitarity bound on the short distance perturbative 
cut-off. Numerical results are provided in Table IVTTT1 for 
Set IV. Typically, we find that results do not depend 
dramatically on the chosen chiral couplings . We take 
r c = 0.5fm as a standard choice. As we see, finite cut-off 
perturbation theory does not work better than our non- 
perturbative results of Sect. E] and in fact requires one 
more counterterm. Actually, this is a perturbative indi- 
cation that NNLO is more important than NLO, casting 
doubts on the convergence of the approach. 

Finally, we have checked that taking the LO to be 
the full OPE potential, Ui„. and the perturbation to be 

{2\ (3) 

V2J 4-U^ as a whole and keeping the cut-off r c > O.lfm 
does not change the results significantly. Actually, the 
perturbative result does not account for the value ob- 
tained non-perturbatively, despite the apparent linear be- 
haviour observed when changing numerically the scaling 
parameter A in the region A > 0.1 (see Fig. I17|) . This 
supports our conclusion that perturbation theory docs 
not compute the slope of As(X) at the origin. In addi- 
tion, even if we disregard the divergence by introducing 
a cut-off, the perturbative calculation does not account 
for the non-perturbative renormalized result. 
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FIG. 19: Dependence of the s-wave function normalization As(r c ), the matter radius r m (r c ), quadrupole moment Qd(r c ) and 
the D-wave probability Po(r c ) in the deuteron in perturbation theory to NLO and NNLO on the short distance cut-off. We 
use Set IV for the NNLO. In any case, the first order solution is fixed to the be normalized and to reproduce the asymptotic 
D/S ratio, rj. We compare with the experimental value and the OPE result. 



TABLE VIII: Comparison of finite cut-off perturbation theory, r c = 0.5fm, with renormalized non-perturbative results for 
the Deuteron properties. We use the non-relativistic relation 7 = ^/2^i np B with B — 2.224575(9). The errors quoted in the 
perturbative calculations reflect the uncertainty in the input parameters 7, r] and As. Similarly, the errors quoted in the TPE 
reflect the uncertainty in the non-potential parameters 7 and r\. We use Set IV of low energy constants ci,C3 and C4. 





r c 


7 (fm- 1 ) 


V 


4s(fin- 1/a ) 


r m (fm) 


Q d (fm 2 ) 


Pd 


(r- 1 ) ( fm- 1 ) 


Non-perturbative 









Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.31(1)% 


0.476(3) 


Uln + U2tt 





Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1)% 


0.447(5) 


Perturbative 




rr(0) 
U 1tt 





Input 


0.02555 


0.8625(2) 


1.9233(5) 


0.2667(1) 


7.14(1)% 


0.484(3) 


u$ + 


0.5 fm 


Input 


Input 


0.862(2) 


1.923(4) 


0.2667(3) 


7.14(1)% 


0.484(3) 


u?J + ug 


0.5 fm 


Input 


Input 


Input 


1.962(2) 


0.2705(2) 


7.12(1)% 


0.484(3) 


Potentials 




Nijmll 




0.231605 


0.02521 


0.8845(8) 


1.9675 


0.2707 


5.635% 


0.4502 


Reid93 




0.231605 


0.02514 


0.8845(8) 


1.9686 


0.2703 


5.699% 


0.4515 


Exp. 




0.231605 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 
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